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ABSTRACT 

Spectral analysis is a powerful tool to investigate stellar properties and it has been 
widely used for decades now. However, the methods considered to perform this kind 
of analysis are mostly based on iteration among a few diagnostic lines to determine 
the stellar parameters. While these methods are often simple and fast, they can lead 
to errors and large uncertainties due to the required assumptions. 

Here we present a method based on Bayesian statistics to find simultaneously 
the best combination of effective temperature, surface gravity, projected rotational 
velocity, and microturbulence velocity, using all the available spectral lines. Different 
tests are discussed to demonstrate the strength of our method, which we apply to 
54 mid-resolution spectra of field and cluster B stars obtained at the Observatoire du 
Mont-Megantic. We compare our results with those found in the literature. Differences 
are seen which are well explained by the different methods used. We conclude that 
the B-star microturbulence velocities are often underestimated. We also confirm the 
trend that B stars in clusters are on average faster rotators than field B stars. 
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1 INTRODUCTION 

Massive stars are bright and short-lived objects. Therefore 
they are relatively easy to observe and remain close to their 
birth place. Their influence on their environment is impor¬ 
tant, from photoionizing HII regions, polluting their environ¬ 
ment with enriched chemical elements, to expelling material 
out of galaxies. But they also represent an important candi¬ 
date, due to their high level of internal radiation pressure, for 
the general study of stellar evolution. Therefore an accurate 
characterization of the fundamental and spectral parameters 
of massive stars is greatly needed in astrophysics. 

One of the most useful tools to investigate massive stars 
is quantitative spectroscopy, which involves the comparison 
of a stellar spectrum with synthetic data created by state-of- 
the-art atmospheric models. There is a wide variety of com¬ 
parison technics presented in the literature. But whether it 
is by fitting spectral li nes “by eye" fSearle et al.||2008[ |, by 
using minimization ( Daflon et al.|20OT l, or by measuring 
equivalent widths i Lefever et al.|2010 l, most of these meth¬ 
ods rely on an iterative process, due to the large number 
of parameters involved in the analysis. For B stars, these 
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parameters are the effective temperature Te//, the surface 
gravity log((?), the projected rotational velocity usin(i), the 
micro- and macro-turbulence velocity ^ and and of course 
the abundances of the chemical elements. The iterative anal¬ 
ysis usually starts by adopting best-estimate values for all 
parameters. It then alternatively refines each parameter in¬ 
dividually while using specific spectral indicators that are 
mostly sensitive to a sub-sample of parameters. For instance, 
in B stars the hydrogen lines H/3 and H7 are mostly sensi¬ 
tive to Tsff and log(gr), while ionized silicon lines are mostly 
sensitive to T^ff, and its abundance. Therefore, one can 
start by assuming initial values for Te//, log(5), and ^ for 
instance, and then use the Balmer lines to determine a new 
log(gr) value. This new result combined with the Si lines is 
then used to infer a new effective temperature and, going 
back to the Balmer lines, to update log(g), and so on until 
all the parameter values are stable. Of course, this does not 
mean that each indicator is completely insensitive to the 
other parameters. On the contrary, there is a highly non¬ 


trivial coupling between all the involved parameters (Nieva 
fc Przybilla||2010 K And due to this coupling, the use of an 


iterative process based on a few spectral indicators can lead 
to a local solution rather than a global solution that takes 


into account all the parameters (Mugnes & Robert 20121. 
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Moreover, depending on the diagnostic lines nsed and on 
the initial parameter values, the same method may not give 
the same results. Also, due to the large number of parame¬ 
ters involved and in order to reduce computation time, many 
works in the literature have focused on a reduced number of 
parameters while assigning fixed values for the others (for 
instance by adopting solar abundance values, or assuming 
^ = 2 km s“^ for main sequence stars and 10 km s“^ for 
evolved stars). Fixed values are often based either on theo¬ 
retical estimates or on results from other works. For exam¬ 
ple, Searle et al.| ( 20081 used the projected rotational veloci¬ 
ties of Howarth et al. 119971. But as discussed in the present 
work, the final results depend on the choice of the param¬ 
eter that is fixed and on the choice of the value adopted. 
Even for a given atmospheric model and a given set of ob¬ 
served spectra, all the simplifications done (i.e. the adoption 
of initial estimated values, the use of fixed parameters, or 
an iterative process with only a few specific indicators for a 
reduced number of parameters) have an impact on the re¬ 
sults and their accuracy and make it difficult to compare or 
to compile results from different works in a consistent way. 
It is therefore important to find a method which does not 
rely on initial estimations of the parameters, that will si¬ 
multaneously constrain all the parameters using most of (if 
not all) the available information contained in an observed 
spectrum, and will work in a reasonable amount of time. 

Thanks to the development of computer power, along 
with its recent success in many scientific fields, Bayesian 
statistic has been more and more used in astrophysics during 
the last decades. Studies found in the literature cover a large 
range of astrophysical topics, from redshift estimations with 
photometry (Benitez 20001 to the study of magnetic fields 


in interstellar clouds (Crutcher et al. 


mological parameter selection (Trotta 


20101 , including cos- 


20071. Among these 


studies, a few authors applied Bayesian statistics to the de¬ 
termination of stellar parameters, hereby illustrating some 
of its advantages. For instance Pont & Eyer (20041 used 
Bayesian statistics to redefine the age-metallicity relation in 
field dwarfs. They showed that most of the scatter found 
in the previously derived age-metallicity relation was due to 
a simplified treatment which, among other issues, did not 
properly take into account systematic biases. [Shkedy et al.| 
(20071 performed a complex and complete study of the im¬ 
pact of both systematic and statistical measurement errors 
on stellar parameters and their uncertainties. More recently. 


Schonrich & Bergemann (20141 developed a method which 


allow the combined use of photometric, spectroscopic, and 
astrometry informations for the determination of stellar pa¬ 
rameters and metallicities. Nevertheless, none of these stud¬ 
ies were applied on massive stars, nor did they performed 
a spectral analysis involving more than three parameters 
at the same time (the spectroscopic part of the work of 
Schonrich & Bergemann 2014 focused only on the effective 
temperature, surface gravity, and global metallicity). 

In this work, we use both a large pre-calculated grid of 
synthetic spectra and small interpolated and tailored grids 
in order to reduce the computation time. Bayesian statistics 
are applied to allow for a logical and numerical connection 
between the two grids as well as to constrain all the param¬ 
eters at the same time while using nearly all the available 
information given by a spectrum. In the present paper, as a 
way to introduce our method and to illustrate the param¬ 


eter interdependency, we focus on four stellar parameters 
(effective temperature, surface gravity, projected rotational 
velocity, and microturbulence velocity). In a following paper, 
we will add individual chemical abundances and the global 
metallicity to the analysis. 

The next section summarises the ingredients and vo¬ 
cabulary used in Bayesian statistics. This is followed by a 
description of our technic and different tests done to demon¬ 
strate its reliability, stellar parameters are obtained for a 
sample of B stars and are compared to other works. 


2 DESCRIPTION OF THE METHOD 


We here briefly describe the key elements of Bayesian statis¬ 
tics involved in our method. For a complete description of 
Bayesian statistics see Gregory (20101. 


2.1 The Bayes Theorem 

2.1.1 Definition 


Let us make an experiment or an observation that allow us 
to retrieve some data D (a stellar spectrum for instance) 
in order to test a given hypothesis H (for example, the hy¬ 
pothesis that a given atmospheric model reproduce correctly 
the stellar spectrum). And let us also assume that we have 
some prior knowledge or information I (the stellar param¬ 
eters given by an approximate spectral type for instance). 
Then by using the Bayes theorem we can assign a plausi¬ 
bility degree (or probability) to the hypothesis H. Namely, 
the Bayes theorem states that the probability P of H being 
true given D and I is called the posterior probability and is 
given by: 


P{H\D,I) 


P{H\I)*P{D\H,I) 

P{D\I) 


( 1 ) 


where P{H\I) is the prior probability which is obtained 
over a previous set of data or from previous knowledge, 
P{D\H, 1) is the likelihood, i.e. the result of the actual data 
D, and where the global likelihood P{D\I) is a normaliza¬ 
tion factor. Simply put, according to the Bayes theorem, 
our final state of knowledge (the posterior probability) is 
the product of what we already knew (the prior probability) 
with what the experiment (the likelihood) tells us. 


2.1.2 Application to a Stellar Spectrum 

In this work we adapt the Bayes theorem to stellar spec¬ 
troscopy considering three main concepts : 

1) Let us consider an observed spectrum of a given star that 
contains several spectral lines. We then artificially divide 
the whole spectrum into smaller spectra, each centered on 
an individual spectral line (Fig. &)• Each smaller spectrum 
can be seen as an individual observation. As the original 
spectrum can be described as a data set D composed of all 
the spectral resolution elements (or flux points) contained in 
the whole spectrum, then each smaller spectrum k can also 
be described as a data set Dk composed of all the spectral 
resolution elements contained in the spectrum k (Fig. 0^)- 
Thus for a given star, we have K independent spectra, where 
each spectrum k is centered on the spectral line k and is 
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Figure 1. Illustration of the data sets, (a): Superposition of an observed spectrum D (black) and synthethic spectrum Mi (dashed red). 
Lighter parts of the spectra are not considered in the analysis, (b): Close up of one sub-spectrum Dj. where each observed and synthetic 
spectral resolution elements are labeled ci„ * and m„ j. i, respectively, with n (1 ^ n ^ Nf^) the number of resolution elements in D^., k 
(1 ^ A: Si K) the number of sub-spectra, and f (1 ^ i Nb^adal) the number of synthetic spectra, (c): Errors e„ i, assuming that the 
observed spectrum is defined as j. ^ -|- 


described by a data set Dk (where 1 ^ k ^ K, and K is 
the total number of lines contained in the original observed 
spectrum). 

2 ) As we want to fit an observed spectrum over a grid of 
synthetic spectra, we can replace the hypothesis H hy a. set 
of hypothesis Mi (with i = 1, 2,..., Nbmodei) representing the 
propositions stating that each synthetic spectrum calculated 
with an atmosphere model M with a set of parameters i is 
true (i.e. correctly reproduces the observation). Then for 
a given data set we can assign a probability for each 
synthetic spectrum. Since the propositions Mi are mutnally 
exclusive, we can write : Yli P[tvIi\D, I) = 1. 

3) Finally each line k (or experiment k) has an associated 
prior information Ik and a posterior probability reflecting 
our state of knowledge respectively before and after doing 
the experiment k. 


We can then rewrite the Bayes theorem for each line k 
as : 


P{Mi\Dk,Ik) 


P{M,\Ik)*P{Dk\Mi,h) 

PiDklh) 


( 2 ) 


If we use the posterior probability given by the line k — 1 
as the prior probability of the line k, then we can define the 
prior information of the line k with Ik = Dk-i, Ik-i (i.e. the 
combined information given by the previous data set Dk-i, 
or line k — 1, and its associated prior information Ik-i)- The 
Bayes theorem then becomes: 


P{M,\Dk,Ik) 


P{Mi\Dk-i, Ik-i) * P{Dk\Mi, Ik) 
P{Dk\Ik) 


( 3 ) 


which can be understood as: the probability that the syn¬ 
thetic spectrum Mi reproduce well the observed line Dk is 
proportional to the quality of the fit (likelihood) between the 
observed line Dk and the corresponding synthetic line from 
the spectrum Mi times the probability that the same syn¬ 
thetic spectrum Mi reproduce the previous line Dk-i with 
the previous corresponding synthetic line. 


Here P{Mi\Dk, Ik) is the posterior probability for one 
synthetic spectrum Mi. When we consider all the synthetic 
spectra at the same time, {P{Mi\Dk, Ik)}i=i..Nb^„^^i^ is 
thus called the posterior probability distribution (the same 
obviously applies for the likelihood and the prior). 

In other words, if we apply this theorem to all the lines 
in the observed spectrum (each time using the last posterior 
probability distribution as a prior for the next line) and con¬ 
sidering all the synthetic spectra, we obtain, in the end, a 
final posterior probability distribution that simultaneously 
takes into account all the constraints given by all the lines 
in the observed spectrum for all the considered parameters. 
Figure shows an example of the process when we apply 
our method on 4 spectral lines using 25 synthetic spectra. 
Note here that the order of the lines to which we apply the 
theorem does not matter. Indeed, if we neglect the normal¬ 
isation factors, the final posterior probability is simply the 
product of the first prior probability with the likelihoods of 
all the lines involved. 


2.1.3 Constructing the Likelihood 

In order to construct a proper likelihood, we need to make 
a few basic assumptions. First, each data set Dk (i.e. line k) 
is composed of Nk independent spectral resolution elements 
Dk = {dn.k} = di,fe,d 2 ,fe, (Fig. [^). Then we can 

do the same for the model spectrum Mi^k = {'fTin,k,i} = 
•••5 (Fi§. &)• Finally let us suppose 

that each observed spectral resolution element can be as¬ 
sociated with a spectral resolution element from a model, 
given a certain flux error d^.k = + e„.k,i (Fig. [^). If 

we suppose the errors {en,k,i} to be independent and repre¬ 
sented by a Gaussian distribution centered on zero with a 
standard deviation of (Jn.k, then we can write the likelihood 
as: 
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P{Du\Mi,h) = ll 


exp 


(^n,fc ^n,fc,z) 


-E 

n 

= Ck exp 




xlk 


( 4 ) 


Finally, when we apply the theorem (Eq. for the 
first time, we assume that all models Mi are equiprobable: 
P{Mi\Do,I(i) = 1/Nbmodei (i-o. we consider a fiat prior dis¬ 
tribution as in Fig. and Do represents no data). This 
means that, at first, we do not give preference to any model 
(or set of parameters). Note here that we are not performing 
a model comparison but a parameter estimation and as we 
use large data samples, the results are not affected by the 
choice of the prior | Gregory ||2010 1. 


The errors can represent the noise in the data but also 
the incompleteness of the model (due to simplifications, ap¬ 
proximations, or unknowed variables) or even some errors 
in the reduction process of the data (bad normalisation for 
instance). To take these sources into account, we define the 
standard deviation of the gaussian distribution as follow: 


(^n,k 


= 'A 


^noise)n,k T i^rnodel)n,k ^ 


( 5 ) 


where {anoise)n,k refers to the errors coming from the noise 
in the data, and {crmodei)n,k refers to the other sources of 
errors (model incompleteness and reduction errors). In this 
work, for simplicity, we assume {anoise)n,k = o-noise, a con¬ 
stant for all the spectrum, defined as the standard deviation 
in regions of the normalised continuum of the observed spec¬ 
trum. And for each line, {amodei)n,k is given by: 


(rr model ')'n,k 



(rr model ) k • 


( 6 ) 


In other words, for each line k, {(Jmodei)k is the minimal 
RMS between the observed line and all the corresponding 
model spectrum Mi^k- With these definitions, an,k is the 
same for all the spectral resolution elements n of a line, but 
is different for each line. Thus: 

rTn.fc — ^noise T {f^modelR ~ (^) 


Note that {o'x)k is a very important variable as it ulti¬ 
mately links physical imperfections and theoretical incom¬ 
pleteness to the uncertainty of the final parameters. And 
with this decomposition, if we were to work with perfect 
data {(Jnoise I —> 0), then the uncertainty of the final param¬ 
eters would be related to the limitations of the atmospheric 
or atomic models used. Inversely, if we have a perfect model 
{{o'modei)k ' —0), the final uncertainties would then be di¬ 
rectly related to the quality of the data. 


2 . 1.4 Global Likelihood and First Prior 

Now, since P{Mi\Dk, Ik) = 1, the global likelihood 
P{Dk\Ik) is then given by: 

P{Dk\Ik) = Y^{P{M,\Dk-i,Ik-i)*P{Dk\Mi,Ik)). ( 8 ) 


And the Bayes theorem (Eq. becomes: 

P{Mi\Dk-i,Ik-i) * Gfeexp 

P{M,\Dk,Ik) = - 


2 


( PiMi\Dk-i,Ik-i) * Gfeexp 


\i,k 


( 9 ) 


2.2 Parameter Estimation 

In this work we focus on the stellar parameters Te// the ef¬ 
fective temperature, log(p) the surface gravity, v sin(i) the 
projected rotational velocity, and ^ the microturbulence ve¬ 
locity. We thus need to retrieve, from the final posterior 
probability mentioned above, an estimation for each param¬ 
eter by considering all the other parameters. For this, we 
use another major feature of the Bayesian statistics, the 
marginalization. This is done by defining the model Mi as: 
the synthetic spectrum computed with the model M using the 
set of parameters Rffi, log{g)i, v sin(i)i, and ^i, which we 
write Mi = Ti,li,Vi,^i, M . Since we suppose our model M 
to be true, our posterior probability thus becomes: 

P{Mi\Dk,Ik) ^ PiP,li,Vi,^i\M,Dk,Ik), (10) 

which can be read as the joint probability for Ti, U, Vi, and 
^i to be true if the model M, the data Dk, and our prior 
knowledge Ik are true. 

Since our parameters have discrete values in our anal¬ 
ysis, we can retrieve the marginal posterior probability dis¬ 
tribution for each parameter by using the marginalization 
equation. Here we give an example for T : 

P{Ti\M, Dk,Ik) = EEE P{T„la,Vl„^e\M, Dk,Ik), 

a h c 

( 11 ) 

which is the marginal probability distribution of Ti con¬ 
sidering all the values of the other parameters. As shown 
in Figure from this distribution we obtain the most 
probable value and the mean value, as well as the asso¬ 
ciated uncertainties which are defined as the values be¬ 
tween which the sum of the marginal posterior probabil¬ 
ity distribution is equal to or greater than a certain confi¬ 
dence level C = P{Ti\M, Dk, Ik). In this work we use 
C = 0.99 which can be seen as a Scr uncertainty (recall that 
P{Ti\M, Dk, Ik) = 1). Note that these uncertainties are 
related to the variations of the other parameters but also, 
as stated in §|2.1.3[ to possible data or model imperfections. 


2.3 Preventing Numerical Biases 

Due to the fact that we want to treat all available lines 
equally, we need to calculate the for the same number of 
spectral resolution elements in each line. Otherwise, broad 
features (such as hydrogen and helium lines in B stars) will 
dominate the solution over the narrow features (metal lines 
for instance) in the posterior probability. 

Indeed, given the same differences dn,k—mn,k,i between 
the observed and synthetic line and the same (j^ for both a 
broad and a narrow line, only the number of spectral resolu¬ 
tion elements will have an impact on the value. The more 
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Data set Dj^ 




/ j 4826 4827 4828 4829 4830 4831 4832 
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Likelihood P(D^|Mj,IJ X Prior P(Mj|Dj^.^,Ii^.i) = Posterior P(Mi|Dt.,y 
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Figure 2. Illustration of the Bayesian analysis step by step. For this example, we consider 4 data sets (or lines) and 25 synthetic 
spectra calculated with only two free parameters: T^ff and The black plus signs in columns 2, 3, and 4 indicate the value used for 
these two parameters. The colour of the surface contours is scaled to the relative probability of each distribution: the dark red regions 
being the most probable and the white ones the less probable. The yellow squares indicate the global maximum of each distribution. 
The first column shows every data set considered, and the second column the likelihood obtained for each data set as in Equation]^ 
The third column represents the prior probability distribution known at each step /c, and the last column shows the posterior probability 
distribution which is the normalised product of the prior and the likelihood of a given step k. The method goes as follow: considering the 
first line k = 1 (in (a)), we calculate its likelihood using Equation]^ (in (b)), and use, as a first prior, a flat distribution (in (c)) assuming 
we have no knowledge of the best pair ^]. We then apply the Bayes theorem (Eq. and obtain the first posterior probability 

distribution (in (d)). For the next data set (fc = 2), we perform the same operation but we use the posterior distribution of the previous 
step (in (d)) as a new prior distribution (in (g)) since our “state of knowledge” has been changed after the first step. In the end, the final 
posterior probability distribution (in (p)) represents the solution combining all the information given by the 4 lines. Note that since the 
method can be summarized as a series of products, the order in which we consider the data sets have no influence on the shape of the 
final posterior distribution (see Fig. 1^. 


spectral resolution elements there are, the higher is the value 
of the x^- Therefore, the x^ of a- broad line will rapidly in¬ 
crease and hence its likelihood will also rapidly tend to zero 
as the models differ from the solution, whereas a narrow line 
will see its likelihood slowly reach zero. Consequently when 
multiplying the prior probability density given by the broad 
line with the likelihood of the narrow line, only the broad 
line solution survives. To prevent this effect, we therefore 
interpolate each line with the same number of flux points 
given by the broadest line in the sample. This interpolation 
will obviously oversample the narrow lines, but it is a good 
way to consider each line equally without losing information 
(since interpolating each line with the number of flux points 
given by the narrowest line would imply a loss of detail in 
the broader lines). The oversampling has also the advantage 
of giving a more realistic weight to the narrow and weaker 
lines; the parameter space for the solution of a weaker line 
is then better defined. 


2.4 Synthetic Spectra 


We use atmospheric models from the metal line-blanketed, 


NLTE, plane-parallel, and hydrostatic code TLUSTY (Lanz 


& Hubeny 20071. With these models, we create synthetic 


spectra by running the program SYNSPEC. More precisely, 
we use the SYNPLOT package generously offered to us by 
Hubeny and Lanz, which contains a TLUSTY pre-calculated 
B-star grid {BSTAR2006) and the latest version of SYN¬ 
SPEC and SYNPLOT which allows the user to interpolate 
more specihc models from the closest BSTAR2006 mod¬ 
els. The BSTAR2006 grid considers 16 values of the ef¬ 
fective temperature, between 15000 K and 30000 K with 
a 1000 K step; 13 surface gravities, from 1.75 to 4.75 in 
log( 5 r) with a step of 0.25 dex; and 6 chemical compositions, 
ZjZQ = 2,1,1/2,1/5, and 1/10 (where Zq is the solar com¬ 
position); a solar helium abundance, He/H = 0.1 by number, 
and a microturbulence velocity of 2 and 10 km s“^. 

For this paper we consider only the solar chemical com¬ 


position, as established by Asplund et al. (20091, and we 
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Figure 3. Exemples of the marginal probability distribution for 
Teff and These distributions (in (a) and (c)) are derived from 
the final posterior probability distribution (in (b)) taken from 
Figure]^. The coloured regions in (a) and (c) indicate the 68%, 
95%, and 99% confidence intervals which can be seen as a Icr, 
2cr, and 3cr uncertainties for the best value of each parameter. 
The dashed red lines show the position of the mean value of each 
parameter defined as {X) = * P(Xj|M, 7^), where X 

is one of the considered parameters. Since the parameter space 
is discrete, the most probable value found for a parameter might 
not be the true value but the closest value used in the calculation 
of the spectra. Thus the mean value can be a better estimation 
of the real value of a parameter. 


use SYNSPEC to create a grid of synthetic spectra (here¬ 
after called BaseGrid) with the same values of the effective 
temperature and surface gravity as those of the BSTAR2006 
grid with an extended range of 7 microturbulence velocities, 
from 0 to 30 km s“^ with a step of 5 km s“^. We then con¬ 
volve these spectra to take into account rotational broaden¬ 
ing with 21 different values of v sin(i), from 0 to 400 km s“^ 
with a 20 km s“^ step, to create our starting grid for the 
analysis: BSTART. Note that the spectra of BaseGrid are 
not convolved with instrumental and rotational broadening 
and are used, in this work, as the basis for the creation of a 
new spectrum or new refined grids by directly interpolating 
the non-convolved spectra (because this is faster than inter¬ 
polating the BSTAR2006 models) to the desired values of 
Rff, log(5), and 

For each star, we begin, as previously mentioned, by 
assuming a uniform prior probability over all the BSTART 
models. We then apply our Bayesian analysis over all the 
selected lines in the spectrum of a star using all the mod¬ 
els of the BSTART grid. This first analysis gives us, after 
marginalization, a reduced parameter space by excluding all 
values outside the 99% confidence level. We then construct 
a rehned grid in this new parameter space and we reana¬ 
lyze the lines assuming, once more, a flat prior. The result 
gives an even more restricted parameter space over which 
we reconstruct another refined grid and redo the analysis, 
and so on. The program stops either when the parameter 
space remains unchanged or when the parameter steps used 
to create a new grid are equal to or smaller than a certain 
limit. Here, as a compromise between accuracy and com¬ 


puting time, we set this limit at 100 K in Tg-ff, 0.03 dex in 
log( 5 r), 5 km s“^ in v sin(i), and 1 km s“^ in 


3 TESTING THE METHOD ON ARTIFICIAL 

STELLAR SPECTRA 

As a way to test the efficiency of the method as well as 
its robustness to noise and uncertainties, we perform three 
tests. 

For the first test, we create an artificial observational 
spectrum using synthetic spectra with the following spec¬ 
ifications: Tf,ff = 21350 K, log( 5 ) = 3.80, v sin( 2 ) = 
45 km s“^, C = 6 km s“^, and solar chemical composition. 
The spectrum is also convolved with a Gaussian instrumen¬ 
tal profile with a full width at half maximum (FWHM) of 
2.3 A. Random noise is also added in order to emulate a 
standard deviation, a noise, that is typical of what we see 
in our actual observations (with Unoise — 0.004 in the nor¬ 
malised flux). Then, we compare the artificial spectrum with 
our models using the method described in § |2.4| But, as we 
are in the situation of a perfect model and imperfect data, we 
set {amodei)k to 0 and we substitute Onoise by xanoise where 
X = 1,2,3,... in Equation!^ This way, while not changing 
the actual noise level in the spectrum, we can overestimate 
its value in the calculation by increasing x, and thus, we can 
test the robustness of the method to noise over-estimation. 
This analysis is then repeated 10 times, each time increasing 
a; by 1. The results are summarized in Figure]^ 





Figure 4. Plot of the stellar parameter values for 10 runs with 
increasing cr-^jUnoise- Individual points and error bars are the 
results of the different runs. The green lines show the real param¬ 
eter values used to create the artificial stellar spectrum (T^ff = 
21350 K, log( 5 ') = 3.80, v sin(i) = 45 km s“^, and ^ = 6 km s“^). 
The orange lines show the results (full lines) and the correspond¬ 
ing uncertainties (dashed lines) from the analysis of the artificial 
spectrum without noise added and for (Txl^noise = 1- The orange 
lines often overlap with the green lines. 


While running the analysis on the artificial spectrum 
without noise added (orange lines in Fig. [^, the program 
returns almost the exact parameter values used to create 
this artificial spectrum (green lines in Fig.[^. There is only a 
slight difference of 10 K for the effective temperature which 
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is well under our minimal uncertainty (100 K). The runs 
for different values of a^/anoiae then give parameter values 
which are all very close to the original values used to create 
the artificial spectrum. 

Increasing the ratio a^/anoiae means overestimating the 
noise level in the artificial spectrum, which should logically 
lead to greater parameter uncertainties since acts as a 
tolerance parameter. Here, the parameter uncertainties in¬ 
crease with roughly as (j„oise/2, as shown in Figure]^ 
This means that even when increasing the flux uncertainty 
by a factor of 2 or 3, there will be nearly no incidence on 
the method accuracy. Note also that, up to ianoiae, the ac¬ 
curacy of V sin(i) is largely higher than the adopted reso¬ 
lution (around 77 km s“'^ for a FWHM of 2.3 A). Usually, 
with classical methods, the uncertainty for v sin(i) with a 
value smaller than half the FWHM should be of the order of 
FWHM/2 since instrumental broadening becomes more im¬ 
portant than rotational broadening under this limit. Theses 
effects are a direct consequence of cumulating the informa¬ 
tion from one line to the others and thus having multiple si¬ 
multaneous constraints on each parameter. Note that when 
we use large values, we find the classical uncertainties 
(close to the v sin(i) value when v sin(i) is below the classi¬ 
cal limit) since, with such a huge tolerance parameter, most 
of the constraints coming from the lines disappears. 



Figure 5. Plot of the relative stellar parameter uncertainties 
Area with increasing (JxNnoiae- For each parameter. Ares 
is the length of the error bars from Figure|^and Amin the length 
of the smallest error bar found for the 10 runs. The thin black 
dashed line is a linear fit of the mean of all the parameter uncer¬ 
tainties. 

Since spectral analysis usually tends to be more difficult 
and less accurate for fast rotating stars due to the blend¬ 
ing of multiple lines, we repeat the first test for a higher 
projected rotational velocity, v sin(i) = 145 km s“^. All the 
other parameters are unchanged (we are using the same lines 
as well for this analysis). The evolution of the parameter un¬ 
certainties is shown in Figurej^ Here, the effect seen in Fig¬ 
ures 1^ and for a low V sin(i) value (compared to the instru¬ 
mental FWHM) disappears while the slope of the fit shown 
in Figureremains the same as in Figure]^ Also, the un¬ 
certainties of each parameter are of the same order as those 
obtained for a slow rotator. As this result is the same for 
an even faster rotator (for instance v sin(i) = 300 km s“^), 
we can conclude that a faster rotational velocity does not 
affect the accuracy on the other stellar parameters when us¬ 
ing our method. Again this is a consequence of cumulating 
the information from one line to the others. 



^ ^notse 


Figure 6. Plot of the relative stellar parameter uncertainties 
Area/Amin with increasing (Tx/crnoise for a higher projected ro¬ 
tational velocity, v sin(i) = 145 km Uncertainty and line 
definitions are the same as in Figure]^ 

For the second test, we again use an artificial spec¬ 
trum (with the following parameters: Teff = 21350 K, 
log(5f) = 3.80, V sin(i) = 125 km s“^, ^ = 3 km s“^, and so¬ 
lar chemical composition) but this time we change the noise 
level added to it in order to create different signal-to-noise 
ratios (SNR), ranging from 25 to 350 per resolution element 
with a step of 25. We estimate the SNR using the der^snr 
program of |Stoehr et al.| ( |2008[ ). On this scale, a SNR of 
250 corresponds to a Onoiae of approximately 0.004 in the 
normalized continuum. For each value of SNR, we create 10 
noisy artificial spectra. For the 140 new spectra thus cre¬ 
ated, we apply our method with the original definition of cr^ 
(Eq.0. This way we simulate a real analysis even though 
we still are in the case of a perfect model and imperfect 
data. Of course, in this case is overestimated by roughly 
a factor of 2, since for a given line fe, {amodei)k ~ Unoiae, 
but as explained above, it has nearly no incidence on the 
method accuracy. The final parameters returned by the 140 
runs are shown in Figure and the overall success rate of 
the method is shown in Figurej^ This test allow us to check 
for the robustness of the method against highly noisy data, 
and we can see that for SNR around 100-125 and higher, 
the method returns correct and accurate results (with un¬ 
certainties equal or close to the minimal uncertainties and 
with a mean success rate of 97.5%). For lower SNR ($; 100), 
the parameter uncertainties become naturally larger but the 
real parameter values are still found nearly 84% of the time 
(within the uncertainties). Considering all the SNR values, 
the overall success rate is around 94%. 

The third test is the same as the second one, except 
that the artificial spectrum is created with the atmospheric 
model ATLAS9 ( [Castelli fc Kurucz|2004| , rather than with 
TLUSTY, using the same stellar parameters. This test is 
close to an analysis of real data since we are in the case of 
imperfect data and inexact models. The results are shown 
in Figure and the overall success rate of the method in 
Figure [TOjCompared to Figure]^ we can see from Figure 
that the best parameters returned by the method for the 
spectrum without noise (orange lines) are different than the 
real parameters (green lines) used to create the spectrum. It 
is especially true for the effective temperature (the other 
parameters are within the minimal uncertainties), where 
TeffBeat = 21930 K wMle TeffReai = 21350 K. This dis¬ 
crepancy is to be expected when comparing full non-LTE 
(i.e. non-LTE atmospheric models with non-LTE radiative 
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transfert) spectra with hybrid LTE/non-LTE (i.e. LTE at¬ 
mospheric models with non-LTE radiative transfert) spectra 
since, non-LTE line profiles tend to be slightly stronger and 
broader than those derived from LTE atmospheres, thus giv¬ 
ing a lower temperature when considering a LTE treatment 


( Huang, Gies fc McSwain|2010 Lanz &: Hubeny|2007| ). Here 
the artificial spectrum is created by a hybrid model while 
the grid of synthetic spectra are calculated using a full non- 
LTE treatment. Therefore our method finds a higher effec¬ 
tive temperature than the one used to create the hybrid 
spectrum. Note that this discrepancy is not taken into ac¬ 
count by cjmodfii since this parameter reflects the model inca¬ 
pacity to correctly reproduce the lines (whatever parameter 
values are considered). Here, our model is able to reproduce 
nearly perfectly most of the lines but for a different effective 
temperature than the one used to create the hybrid lines. 
With this in mind, if we consider the best parameter values 
found for the hybrid spectrum without noise as the refer¬ 
ence values (rather than the real values used to create the 
hybrid spectrum), we find a mean success rate of nearly 93% 
(Eig. |10[ ) which is equivalent to the success rate found in the 
previous test (Fig. [^. Note that if we consider the original 
parameter values, the mean success rate for all the parame¬ 
ters is around 75%, and around 96% when we exclude Tiff¬ 
in tests 2 and 3, if we multiply the derived parameter uncer¬ 
tainties by a factor 1.5 or 2, the overall success rate (which 
is around 94% in both case) jumps to nearly 98% and 99%, 
respectively Thus, in terms of classical statistics, our results 
are given with an accuracy close to a 2(j uncertainty. 




SNR 


Figure 8. Success rate of the method, to derive the correct pa¬ 
rameter values (within the uncertainties) used to create the orig¬ 
inal artificial spectrum, as a function of the signal-to-noise ratio. 
Percentages are mean values gathered from all the parameters for 
each SNR value. 



Figure 9. Same test as in Figure[^when using an artificial spec¬ 
trum created by the atmospheric model ATLAS9. The green lines 
show the real parameter values used to create the original artificial 
spectrum = 21350 K, log(g) = 3.80, v sin(i) = 125 km s“^, 

and J = 3 km s“^). The orange lines show the best param¬ 
eters returned for the original spectrum without noise added 
{T^ff = 21930 K, log(g) = 3.81, v sin(i) = 125 km s“^, and 
^ = 4 km s“^). The orange line overlaps the green line in the case 
of V sin(j). 


4 TESTING OUR METHOD ON REAL 
STELLAR SPECTRA 

4.1 Observational Data 


Figure 7. Plot of the stellar parameter values as a function of the 
signal-to-noise ratio. For each SNR value, 10 artificial spectra with 
a random noise structure are tested. The green lines show the real 
parameter values used to create the original artificial spectrum 
[T^ff = 21350 K, log(g) = 3.80, v sin(i) = 125 km s“^, and ^ = 
3 km s”'^). The orange lines show the best parameters returned 
for the artificial spectrum without noise added. The orange lines 
often overlap with the green lines. 


The selection criteria for our sample were as follow: 1) 
only mid-B to early-B type stars, since the pre-calculated 
BSTAR2006 atmospheric models only cover the range of ef¬ 
fective temperature for these subtypes; 2) relatively bright 
stars to allow a high SNR; 3) no binaries or peculiar stars to 
avoid complications while testing our method; 4) relatively 
close objects in order to work with a solar chemical abun¬ 
dance approximation (chemical element abundances will be 
the topic of a following paper); and 5) well known objects 
to allow the comparison of our method with other works. 
The main studies used here for the comparison are from 


Huang, Gies & McSwain (2010 hereafter HGIO), Lefever 
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Figure 10. Success rate of the method, to derive the correct 
parameter values (within the uncertainties), as a function of the 
signal-to-noise ratio. Here the correct parameter values are ei¬ 
ther the values used to create the hybrid artificial spectrum {Real 
param), or the best values found when analysing the artificial 
spectrum without noise {Best param). Percentages are mean val¬ 
ues gathered from all the parameters for each SNR value. 


et al. (20101, Takeda et aL| (20101, 

Searle et al. (20081, 

Markova & Puls ||2008l, Daflon et al. ( 

20071, Huang & Gies 


||2006a|b| hereafte r HGOG), |Lyubimkov et al.| ( |2005| ), and 


Andrievsky et al. ( 1999|). 


We test our method on a sample of 54 mid-B to early-B 
stars, where 38 are from the field and 16 from clusters. The 
field stars are all nearby dwarf and giant B stars, with a visi¬ 
ble magnitude 3.8 ^ 14 ^ 9 . 7 . The cluster stars are members 
of two open clusters: 8 B stars from NGC1960 (M36) with 

8.8 14 5 ; 10.8, and 8 B stars from NGG884 (x Persei) with 

8.9 ^ 14 ^ 11.6. These two clusters have a colour excess 
E{B — 14) = 0.22 and 0.56, respectively, and due to their 
proximity, are assumed to have a metallicity close to solar 
(Kharchenko et al. 2005 I. The spectra have been collected at 
the 1.6 m telescope of the Observatoire du Mont-Megantic 
at 3 epochs: February 2011, August 2013, and December 
2013. They cover the wavelength range from 3500 to 5500 A 
with a spectral resolution of 2.3 A and a mean SNR of 250 
per pixel in the continuum (with very little variation from 
the blue to the red end). We select all the observed lines 
available for the spectral analysis, except when they are not 
reproduced by the models, when they are barely distinct 
from the noise, or when they are superimposed on bad CCD 
pixels. Tablelists all the major useful lines. Note that not 
all these lines are available for every stars since they may 
appear or disappear depending on the spectral type. Note 
also that due to our modest resolution, and sometimes due 
to strong rotational velocities, more minor lines are included 
in the analysis. 


4.2 Results and Discussions 

We apply our method using most of the available lines for all 
54 spectra. Figures [TT] and[T 2 |show 4 examples of the best fits 
obtained with a small and a high rms, respectively. All the 
other quality fits are between these extreme examples. For 
all stars, the fit is globally good when considering the whole 
spectrum while it may vary from excellent to decent when 
considering each line individually. For instance, F1D215191 
in Figure ED has a global rms of 0.0038, while individual 
lines have a rms between 0.0011 and 0.0085. This is because 
the method returns the best solution for all the lines at 


Table 1. Principal spectral lines used for the analysis 


Species 

Wavelength (A) 


H I 

3712, 3722, 3734, 3750, 3771, 3798, 3835, 
3970, 4102, 4340, 4861 

3889, 

He I 

3785, 3820, 3867, 3872, 3927, 3936, 4009, 

4026, 


4121, 4144, 4169, 4388, 4438, 4471, 4713, 
5016, 5048 

4922, 

C H 

4267, 4374, 4619, 4619, 5133, 5145, 5648 


N H 

3995, 4044, 4228, 4237, 4447, 4601, 4607, 
4803, 4994, 5667, 5680 

4631, 

O H 

3912, 3954, 3982, 4070, 4079, 4085, 4153, 

4185, 


4277, 4304, 4317, 4320, 4367, 4415, 4501, 
4642, 4649, 4662, 4676, 4705, 4907, 4943 

4591, 

Mg II 

4481 


Al HI 

4513, 4529, 5697 


Si H 

3856, 3863, 4131, 5056 


Si HI 

4553, 4568, 4575, 4683, 4829, 5740 


Si IV 

4089 


S H 

4294, 4525, 4816, 5032, 5103, 5321, 5433, 5454 

S HI 

4254, 4285, 4362 


Fe H 

5169, 5260 


Fe III 

5074, 5087 



the same time and not the best solution for each line taken 
individually. 

Concerning NGC1960 SAB30 (Fig. 12 upper panel), the 
poor quality of the fit (by far the worse case in our sample) 
is mainly due to the fact that this star has an effective tem¬ 
perature beyond the reach of our model. Its effective tem¬ 
perature is most likely around 13000 K (HG06), while our 
grid can only go down to 15000 K. This is also the case 
for NGC1960 SAB23 (Fig. 12 second panel) but to a lesser 
extent, since its effective temperature seems to be closer to 
15000 K (HG06) than for SAB30. 

For HD40160 (Fig. lower panel), we encountered 
some difficulty during the reduction process, due to the pres¬ 
ence of humidity spots on the CCD at the time of the ob¬ 
servation, making the H/3 line unusable, altering the shape 
of H 7 (i.e. a shallower centroid), and reducing the number 
of usable lines in general. We still performed the analysis 
of this star as an extreme test for the sensitivity of our 
method to the number and quality of the lines used. While 
our results for HD40160 are not extremely different from 
those found in the literature, they present one of the largest 
discrepancy of our results with other works, and therefore 
are not considered in the remaining of this paper. In par¬ 
ticular, we find Te// = 2 OOIOI 240 K, log((/) = d.Oblo^oL 

V sin(i) = 128lg km s“^, and ^ = 5^2 km s“^, while 
HG06 found = 16453 ± 250 K, log( 5 ) = 3.735 ± 0.027, 

V sin(i) = 115 ± 35 km s“^, and ^ ~ 2 km s“^. 

Figure compares our results with those from other 
groups. All the numerical values can be found in Tables [BT] 
(for field stars) and Table B2 (for cluster stars), where the 
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Wavelength [ A] 


Figure 11. Examples of the best fits obtained for four B stars in our sample. These fits have the smallest averaged rms of our sample 
(from top to bottom, the averaged rms are 0.0043, 0.0042, 0.0041, and 0.0038). Observed spectra are in black, best synthetic spectra are 
in red, and the offset rms are in green. Regions of spectra not used for the analysis are in lighter colour. 


uncertainties given by our method (as explained in § 2.21 
are generally smaller than those found in the literature, be¬ 
cause of the cumulative effect of our method (as explained 
in Fig.[^. While most of the time there is a good agreement 
between our results and those from other works, there are 
some important discrepancies. These are discussed in the 
following subsections. 

4-2.1 Chemical Abundances 

In this paper we wants to do a study of the impact of the 
spectral synthesis method used on the stellar parameters, 
therefore it is important that we first demonstrate that the 
discrepancies in Figure are not associated to the chem¬ 
ical abundances considered by the different analysis. Some 
of the other works used here for comparison derived the 
abundances of a few chemical elements at the same time 
that they determined the basic stellar parameters, while we 
are considering a solar abundance, based on |Asplund et ah] 
(20091, as the standard chemical composition for the mod¬ 


els for all the stars in our sample. Nevertheless, the vari¬ 
ations of the abundances found by the other authors are 
small, the overall metallicity is solar, and this distinction is 


not the main reason for the discrepancies found in Figure 
|13| To demonstrate this, we repeat the analysis of the stars 
in our sample using the abundances derived by the other 
works. Specifically, we adopt the abundances of the follow¬ 
ing chemical elements: He and Si from LEFIO (Table [B3| ); 
O and Ne from TAKIO (Table [B4| ); C, N and O from SEA08 
(Table B5l; He and Si from MAR08 (Table B6|; C, N, O, 
Mg, Al, and Si for DAF07 (Table [B7| ; He and Mg from 
LYU05 (Table B8l; and C and N from AND99 (Table [B^ , 
for the stars that we have in common. HGIO and HG06 did 
not derive the abundances for any specific elements but they 
used the atmospheric model ATLAS9 with the solar abun¬ 
dances from Grevesse fc Sauvalj (19981, as we also redo for 
the corresponding stars (Tables B10| and |B11 1 . Note that 
HGIO and HG06 derived Te// and log(^) using ATLAS9, 
but determined v sin(i) using TLUSTY models and as such, 
one should consider only Figure for a proper comparison 
of V sin(i). Figure [I^ shows the comparison of our new re¬ 
sults, with the other papers, obtained when we are using 
their abundance data. When comparing Figures and 
we see a small and marginal difference between the two set 
of results (except for the stars studied by HGIO which we 
discuss in more details in §|4.2.2| We believe that the reason 
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Figure 12. Examples of the worst fits obtained for four B stars in our sample. These fits have the highest averaged rms of our sample 
(from top to bottom, the averaged rms are 0.05, 0.011, 0.0098, and 0.0088). Colours are as in Figure [E 


for the small change between the two figures is mainly re¬ 
lated to the fact that the element abundances derived by the 
other authors were tailored to a few lines only and have lit¬ 
tle impact on the other parameters obtained simultaneously 
by fitting the whole spectrum as is done with our method. 
As discussed in the following subsections, we find that the 
impact of the method used and the choice of the diagnos¬ 
tic lines is more important to explain the discrepancies seen 
in Figures or A complete analysis with our method, 
including abundances determination and their effect on the 
basic stellar parameters, will be done in a following paper. 


4-2.2 Effective Temperature and Surface Gravity 

Figure |13| or |14| often shows a large scattering for the ef¬ 
fective temperature and the surface gravity. Among others, 
we see a quasi-systematic difference between our estimates 
of Te// and log( 5 ) and those of HGIO and HG06. We be¬ 
lieve that the main reason for this is related to the fact that 
these authors used only H 7 as a diagnostic line. For the 30 
stars that we have in common with HGIO and HG06, we 
perform our analysis using the atmospheric model ATLAS9 
and the solar chemical composition of|Grevesse fc Sauval| 


(1998 1 with all the available spectral lines (these results are 
presented in Fig. |14[ ), but also, with only H 7 as a diagnos¬ 
tic line. Figure |15| shows for the 30 stars the effect on Te// 
and log((;) of the number of diagnostic lines used, as well 
as the effect of the different models. From the upper panels 
of Figure |15[ we see that a change of atmospheric model 
or chemical composition (TLUSTY with the solar composi¬ 
tion of |Asplund et aL]|2009|in red versu s ATLAS9 and the 
solar composition of Grevesse fc Sauval||1998 in green) has 
little impact on the discrepancies between our results and 
those of HGIO and HG06: the discrepancies remain impor¬ 
tant. Whereas the lower panels shows that considering only 
H 7 as a diagnostic line drastically reduces the variations in 
Te//, and greatly affect those in log(gr). For instance, the 3 
stars in Figure [T3| with a large difference in Te//, when com¬ 
pared to HGIO, are the most extreme examples of this effect 
(they are the stars number 23, 24, and 25 in Figure [T^ . It is 
clear in this case that the number and choice of diagnostic 
lines used has a greater impact on the determination of Te// 
and log((/) than the choice of the atmospheric model or the 
solar chemical composition. 

Furthermore, the choice of the diagnostic line, when 
only one is used, also affects the results. Figure shows 
the best values of Te// and log(g) that we obtain for differ- 
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Log{g) [cgs] (literature) 



vsin{i) [km s'^] (literature) 



Figure 13. Comparison of our results for log(g), v sin(i), and ^ with other works from the literature. Atmospheric models from 

TLUSTY with the solar composition from Asplund et al (1999) have been used. Open diamonds are cluster stars and the other symbols 
are field stars. Mean error bars for these data are shown in black. The black full line indicates a one-to-one match. Note that some of 
our stars have been studied in more than one paper and thus appear multiple times in the graphics. The abbreviations in the legend 
correspond to the following papers: H GIO for|Huang, Gies Sz McSwain| | |2010ll; LEF10|Lefever et alT] ) |2010[l; TAK10|Takeda et al.| l |2010[ |; 
SEA08 [Searle et al.| l|2008| |; MAR08 [Markova &: Puls| l|2008||; D ATO7 [Daflon et al.| ) |2007| l; HG06 [Huang fc Gies| | |2006a|b^ ; LYU05 
[Lyubimkov et al.j (|2005|l;~^d AND99 for [Andrievsky et al.|||19^. 


ent values of v sin(i) when only H/3, Hy, or H5 is considered 
for the analysis. This analysis is done using an artificial spec¬ 
tra with Te// = 21350 K, log(o] = 3.8, ^ = 3 km s“^, and 
V sin(i) = 15 km s“^. In Figure [l6| we can see that each line 
has a different sensitivity over v sin(i) and that two lines 
rarely give the same results for a given v sin(i). Moreover, 
Te// and log(gr) given by Hy are lower than the real values 
by 2350 K and 0.3 dex, respectively, when v sin(i) is close 
to its real value. These differences are of the order of the 
quasi-systematic difference seen between HGIO, HG06, and 
our work. 


When comparing our results for the evolved stars with 
those obtained by Searle et al. (20081 in Figure [l^ we find 
similar temperatures for the stars below 20000 K, and less 
scatter than [Searle et al.'| ( |2008| ) for the stars above (the effec¬ 
tive temperature goes from 23500 K to roughly 25000 K in 
our work, while the authors found temperatures going from 
21000 K to 30000 K). The log{g) values are, within the un¬ 
certainties, quite identical, with log( 5 r) below 3 in our case. 


We also find higher microturbulence velocities and a larger 
scatter for the rotational velocity. Note that performing our 
analysis with the C, N, and O abundances derived by[Searl^ 


et al. (20081 has little to no impact on the observed discrep¬ 


ancies (as shown when comparing Figures 13 and 141. The 


reason for the differences between our results and those of 


Searle et al. (20081 may be related to two important effects: 


1) First of all, Searle et al. (20081 used a combination of 


TLUSTY and CMFGEN in order to include wind effects 
which we are not doing. For instance, in B supergiant, the 
Balmer lines could suffer from wind contamination, where 
photons emitted by hydrogen atoms in the wind “fill in” the 
absorption profiles thus making them shallower than any 


profiles calculated with a photospheric model (Searle et al. 
[2008p . It is especially true for Ha and H/3 and to a lesser ex¬ 
tent for Hy and the other Balmer lines. Other photospheric 
lines can also be influenced by winds, such as the Si IV A4089 
line and the Si III AA4552,4568,4575 multiplet ( [Hillier et al.[ 
2003 Dufton et al.|[2005 Searle et al.||2008 K But since this 
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Log(g) [cgs] (literature) 
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vsin(i) [km s'^] (literature) 


I [km s '] (literature) 


Figure 14. Comparison as in Figure |13| but with the individual chemical abundances fixed, in our analysis, to the values given by 
other papers. The modified chemical elements are as follow: He and Si for a comparison with LEFIO; O and Ne for TAKIO; C, N, and 
O for SEA08; He and Si for MAR08; C, N, O, Mg, Al, and Si for DAF07; He and Mg for LYU05; and C and N for a comparison with 
AND99. No individual chemical abundance modifications were considered for the stars in common with HGIO and HG06, but we used 
the atmospheric model ATLAS9 with the solar abundances of |Grevesse &: Sauval] l |1998[ l in this case. Note that we find some stars with 
15000 K since the grid calculated with ATLAS9 has a broader range in T^ff than our grid created with TLUSTY. Also, while 
HGIO derived Te// ^ud log(g) using ATLAS9, these authors used TLUSTY to derive v sin(i), the comparison for this parameter is thus 
irrelevant here. 


influence seems to be model dependant, it is therefore dif¬ 
ficult to correctly quantify it. Nevertheless, it is not really 
surprising that among the 8 program stars that we have in 
common with jSearle et al.] ( |2008| ), the 3 stars that exhibit the 
strongest differences in Te// and log((?) are the ones with the 
strongest mass-loss rate (according to Searle et al. 2008 mea¬ 


surements). These 3 stars are HD192660, HD213087, and 
HD190066. Respectively, their mass-loss rates are 5, 0.7, and 
0.7 X 10“®Mq yr~^, their effective temperature differences 
are 4290, 3250, and 2370 K, and their log(( 7 ) differences are 
0.3, 0.27, and 0.07 dex. 


2) Furthermore, Searle et al. (20081 used, in an iterative way, 


different diagnostic lines for different parameters each time, 
while fixing the other parameters and adopting a rotational 
velocity found in the literature. While vastly used, due to 
its simplicity and swiftness, this iterative method may pro¬ 
vide a local maximum in the probability space rather than a 
global maximum because it considers each parameter inde¬ 


pendently. For instance, [Searle et al.| | |200^ used Si IV A4089 
and Sim AA4552,4568,4575 as primary diagnostic lines for 
the effective temperature with fixed values of v sin(i) and 
Then they used the effective temperature obtained and 
the Hq line to determine log)^). And finally they adopted 
their pair [Te//, log(g)], and a value of v sin(i) found in 
the literature, to determine a new microturbulence veloc¬ 
ity using the Si III AA4552,4568,4575 lines. The whole pro¬ 
cess was repeated until small changes in the parameters 
were seen. With this iterative method, the effective tem¬ 
perature returned by the silicon lines greatly depends on 
the values used for the microturbulence and rotational ve¬ 
locities. And log( 5 ) obtained from Hy also depends on the 
effective temperature. Finally, the microturbulence is signif¬ 
icantly affected by Te// and \og{g). To illustrate the impact 
of an iterative method, we create another artificial stellar 
spectrum with the following specifications: Te// = 25000 K, 
log(gr) = 2.8, V sin(i) = 100 km s“^, and ^ = 14 km s“'. 
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Figure 15. T^ f f and log (g) variations for different models and different number of lines for the 30 stars in common with the sample of 
|Huang, Gies McSwain|(|2010|l and|Huang Gies|(|2006a||. The variations are defined as AT^ff = T^ff (this work) —T^ff {Huang et al.) 
and A\og{g) = \og{g){this work) — \og{g){Huang et al.). Upper panels: Variations when considering all the lines but different atmospheric 
models and solar chemical composition. Lower panels: Variations when considering all the available lines or only H 7 . 



Figure 16. The most probable pair [T^ff, log(^)] as a function 
of V sin( 2 ) for different lines (as indicated in the plot) for the 
artificial spectrum. The parameters used to create the artificial 
spectrum are: Tg// = 21350 K, log(p) = 3.8, ^ = 3 km s~^, and 

V sin( 2 ) = 15 km s“^. The behaviour as a function of v sin('i) 
is shown projected in the plane [Teff, log(p)]. The likelihoods 
are marginalized over all ^ values. Note in particular the most 
probable pair [T^ff 19000 K, log(^) = 3.5] given by H 7 when 

V sin( 2 ) ^ 150 km s~^. 


We then calculate the posterior probability for the 4 silicon 
lines from this artificial spectrum with the BSTART grid. 
After the marginalization over the log( 5 f) domain, we find 
the most probable effective temperature for each given value 
of ^ and V sin( 2 ) (Fig. 17 left panel). We then perform the 
same operation for H 7 and marginalize over v sin(i) and ^ 
(since both parameters have little influence on the returned 
log( 5 f) value, showing the marginalized probability is enough 
in this case) to obtain the most probable log( 5 ) estimate 
given each value of Teff (Fig. |17[ middle panel). And fi¬ 
nally we redo the same analysis, but this time for only the 
3 lines Si III AA4552,4568,4575 and with a marginalization 
over V sin(i) in order to find the most probable value of 
^ for each value of Teff and log(g) (Fig. 17 right panel). 
To conclude, we can see that the chosen silicon lines al¬ 
most always overestimate the effective temperature (except 
when ^ = 10 km s“^ and v sin(i) ^ 200 km s“^). Figure 17 


also shows a clear anti-correlation: the effective temperature 
decreases as the rotational velocity increases. Thus if we 
were to underestimate (respectively overestimate) v sin(i), 
we would overestimate (respectively underestimate) Teff. 
As log(g) increases with the effective temperature, we would 
therefore overestimate (respectively underestimate) log(g). 
Finally, we would overestimate or underestimate the micro¬ 
turbulence velocity depending on how we estimate log((/). 
This behaviour is what we observe between our results 
and those of Searle et al. (2008 1 for the stars hotter than 
18000 K. For the cooler stars, Searle et al. (20081 used 
Si III AA4128,4130 rather than SiIVA4089, which seems to 
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Figure 17. Illustration of the interdependency between T^ff, log(g), and ^ when individual diagnostic lines are considered in an iterative 
method. The parameters of the artificial spectrum used in this exemple are: T^^ff = 25000 K, log(g) = 2.8, v sin(i) = 100 km s“^, and 
^ = 14 km s“^. Left Panel-. Most probable effective temperature as a function of v sin(i) for 3 different microturbulence velocities 
using Si IV A4089 and Si III AA4552,4568,4575. Orange squares and thin orange dashed lines indicate the extreme solutions found when 
underestimating and overestimating v sin(i) compared to the parameters used to create the artificial spectrum (thin green dotted lines 
in all panels). Middle panel: Most probable log((;) values as a function of after marginalizing over v sin(f) and g while using Hq. 

Orange squares and orange dashed lines indicate the solution found when overestimating and underestimating T^^ff from the left panel. 
Right Panel: Most probable microturbulence velocities as a function of T^ff, after marginalizing over v sin(i), for 3 different log(g) 
estimates. The full green, dashed blue, and dot-dashed red lines represent the possible values of log(g) found in the middle panel when 
overestimating and underestimating T^ff. The thin orange dashed lines indicate overestimated and underestimated f from the left 
panel. 


induce a smaller effective temperature discrepancy with our 
method. 

Results from the other authors in our comparison list 
are in a reasonably good agreement with our results even 
though a wide variety of technics and methods were used. 
Takeda et al. ( 2010[), [Lyubimkov et S? (20051, and An¬ 


drievsky et al. ( 1999|) used the uvbyj photometric system 


to derive Te// and log((/). Our results are quite consistent 
with theirs. Daflon et al. ( |2007 1 derived Te// using the UVB 
photometric system and obtained log( 5 ) by adjusting the ob¬ 
served Hy line wings with models from ATLAS9. Note that 
we do not observe the same discrepancies in log((?) with 


Daflon et al. (20071, as seen between our results and those 


of HGIO, since on the one hand, Daflon et al. (20071 used 
Tf.ff values close to ours, and on the other hand, they did 
not fit the core of the Hy line where non-LTE effects are 
stronger. Markova & Puls (20081 used a similar approach 
than [Searle et al.| ( |2008| ): models including stellar wind, 
log(gr) from Balmer lines, and Te// from Si lines. Our results 
agree well with theirs since the stars we have in common 
have weak stellar winds and Te// below or around 18000 K. 
Finally, Lefever et al. (20101 used a semi-automated iter¬ 
ative method which is close to our simultaneous approach 
because groups of parameters are derived using a large num¬ 
ber of lines simultaneously. Our respective results are in a 
reasonable agreement considering that [Lefever et ah (20101 
used models including stellar wind (FASTWIND) and an 
iterative method. 


4-2.3 Projected Rotational Velocity 

A substantial difference in Figure |13[ between our results 
and those published in the literature, concerns v sin(i). The 
most important discrepancy is found for 4 cluster stars that 
are also in the sample of HGIO (as mentioned earlier, they 
also used TLUSTY to derive v sin(i)). We indeed derive 
higher v sin(i) values (around 300 km s“^) by a factor of 



Wavelength [A] 



Figure 18. Comparison of the best fits for the lines HelA4471 
and Mgll A4481 for 2 stars with large rotational velocities. Indi¬ 
vidual data points with error bars represent our observations. In 
full red: models based on our best set of parameters. In dashed 
green: models using the parameters of HGIO. The v sin( 2 ) values 
are: 288 km (red) and 77 km s“^ (green) for BD-l-34 1103s; 
and 315 km s“^ (red) and 199 km s“l (green) for SA058329. The 
respective complete sets of parameters are listed in Table [B^ 
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Wavelength [A] Wavelength [A] 



Figure 19. Comparison of the fits for 4 helium lines obtained with our best set of parameters and those from [Lyubimkov et al.| l |2005[ l, 
for HD184171. Individual data points with error bars represent our observations. In full red: our best mo del with T^ff = 16890 K , 
log(g) = 3.89, § = 3 km s"'^, and v sin(i) = 3 km s“^. In dotted green: model using the parameters of [Lyubimkov et al.| | |2005| : 
Tf,ff = 16100 K, log(g) = 3.62, ^ = 0 km s“^, and v sin(i) = 29 km s“^. 


at least 1.5 for BD+34 1103s, NGC1960 SAB20, NGC1960 
SAB30, and SA058239. Of these stars, SAB20 and SAB30 
have effective temperatures below our limit of 15000 K and 
the lowest SNR (around 200 and 75 respectively) objects of 
our sample. Consequently, the best fits for these stars are 
among the worst in our sample (Fig. 121 even when using 
ATLAS9 models, we thus consider that the associated re¬ 
sults are not reliable. On the contrary, the other two stars, 
BD-l-34 1103s and SA058239, have excellent SNR and rms 
(see Fig. 11 for SA058239). Here, the differences in v sin(i) 
cannot be explained simply either in terms of the number 
or the choice of lines, or in terms of the atmospheric mod¬ 
els used, since HGIO also used TLUSTY and SYNSPEC to 
determine v sin(i) and we could not reproduce their results 
using the same lines. One plausible explanation would be 
the difference in the treatment of v sin(i) (HGIO used a grid 
of models based on a three-dimensional parameterization of 
Tsff, log(g), and the cosine of the angle between the sur¬ 
face normal and the line of sight) but, in this case, there 
should be a systematic effect between our derived values 
and theirs for all the cluster stars that we have in common, 
which we do not observe. Moreover, when we reproduce the 
lines HelA4481 and MgIlA4481, used by HGIO to derive 
V sin(i), we clearly obtain an excellent fit of our own data 
with our parameters (Fig. 181. We are thus confident of our 
derived parameters for these two stars even though we have 
no clear explaination for these discrepancies. 


Interestingly, while the differences between our results 
for V sin(i) and those of HGIO and HG06 are generally 
greater for cluster stars than for field stars, and even though 
our sample is not large, we still find the same conclusion as 
HG06, i.e. cluster stars are globally faster rotators than field 
stars. From our results, cluster stars are even faster rotators. 


Another issue concerning v sin(i) in Figure 13 


14 


appears at relatively low rotational velocities (v sin(i) ^ 
100 km s“'^). In this case, there is a rather large scatter 
between our results and those from other authors. Some 
of the largest differences can either be explained by our 


limite d range in temperature (for example, Takeda et al. 
2010 found V sin(i) = 15 km s“^ for HD185330, while we 


find V sin(i) = 115 km s“ because the effective tempera¬ 
ture is most certainly below 15000 K), or by the parameters 
and spectral lines considered. For instance, the v sin(i) val¬ 
ues of [Lyubimkov et al.j ( |2005| ) are mostly between 15 and 
50 km s“^, while ours are between 0 and 10 km s~^ with 
an uncertainty around 10 km s“^ for the same stars. Lyu- 


[bimkov et al.[ ( [2005[ ) derived their velocities using fixed val 
ues of Teff and log( 5 r), free values of v sin(i), and helium 
abundances, and considered 6 helium lines, while we use so¬ 
lar helium abundance, free values of Te//, log((l), v sin(i), 
and and more spectral lines. While helium abundance has 
an obvious impact on the shape of helium lines, relatively 
small changes in Te//, log(ff), and ^ also have an influence 
on the resulting v sin(i) value. For instance, when consid- 
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ering HD184171 with = 16100 K and logg = 3.62, 


Lyubimkov et al. (20051 found a solar helium abundance, 
V sin(i) = 29 km s“^, and ^ = 0 km s“^, while we ob¬ 
tain Tf,ff — 16890 K, logg = 3.89, v sin(i) = 3 km 
and 5 = 3 km for a solar helium abundance. Figure 


19 


shows the models for these 2 sets of parameters for 4 helium 


lines (that are among the 6 helium lines used by Lyubimkov 
et al. ][2005l ). From this figure, it is clear that, with our set 


of parameters, a higher value of v sin(i) would not improve 
the fits. Note also that our parameters are derived from the 
analysis of more than these 4 helium lines (27 in this case), 
and as such are more constrained, thus explaining the small 
uncertainties. Indeed, if we perform an analysis using only 
these 4 helium lines, with Te//, log((;), and 5 fixed to their 
best values, the uncertainty for v sin(i) goes roughly from 
±10 to ±20 thus reducing the discrepancy between our re¬ 
sults and those of Lyubimkov et al. (20051. 

There is also an interesting point concerning low 
rotational velocities and their respective uncertainties. 
Indeed, HGIO noted that depending on the resolution of a 
given spectrum, and for v sin(f) less than half the FWHM 
of the instrumental broadening function, the instrumental 
broadening dominates the rotational broadening, implying 
uncertainties of the order of ±FWHM/2 for low v sin(i) 
values. In the work of HGIO, FWHM « 0.8 and 1.5 A 
resulting in uncertainties of 29 and 48 km s“^, respectively, 
for V sin(i) below these values. Our data have a FWHM of 
2.3 A, and thus the theoretical limit is around 80 km s“^. 
But, as our method cumulates the constraints of a large 
number of lines, it allows us to significantly reduce this 
limit. Figure shows, for each star, the maximum un¬ 
certainty of V sin(i) as a function of v sin(i). Here, we see 
that the method generally and naturally returns larger 
uncertainties for v sin(i) 80 km s“^ than for higher 
values, but with a maximum uncertainty (for low v sin(i)) 
of roughly 30 km s“^ well below the limit of 80 km s“^. 
As a further demonstration, we apply our method on 
higher spectral resolution data of two stars in our sample, 
obtained with a different grating at the Observatoire du 
Mont-Megantic (in February 2011). These two stars are 
HD77770 and HD44700, whose spectra have an effective 
resolution of 0.9 A and a wavelength range spanning 
from 3600 to 4400 A. The results for Te//, log(g), and 
5 obtained for these two stars are the same (within the 
uncertainties) at moderate and high resolution. But for 
V sin(f), HD77770 gives OlJ^ km s“^ at high resolution, 
rather than 14^j4 km s“^ at moderate resolution, and 
HD44700 gives v sin(z) = 4 ^ 4 ° km s“^, at high resolution. 


rather than 28 I 22 km s ^ 


at moderate resolution. 


4-2.4 Microturbulence Velocity 


Finally, one of the most obvious difference in Figure or 
[141 lies with the microturbulence velocity but it is simply 
explained by the fact that, in most works, the microturbu¬ 
lence velocity is not a free parameter. Usually 5 is fixed at 
about 3 km s“^ and 10 km s“^ for main-sequence stars and 
for evolved stars, respectively (HGIO; [Takeda et al. 20101. 
And, when the microturbulence velocity is not fixed, it is 
estimated with only one or two lines in the spectrum ||Lyu- 


bimkov et al.|[2005l or with only one species (Lefever et al. 



Figure 20. Maximum uncertainty for v sin(i) as a function of 
V sin(i) for our whole sample. 


|2010p . When considering all the available lines to fit the 
microturbulence velocity simultaneously with the other pa¬ 
rameters, we find a higher estimate than what is gener¬ 
ally used. We obtain a mean microturbulence velocity of 
6.1 km s“^ with a standard deviation of 2.6 km s“^ for the 
main-sequence stars, and a mean value of 21.4 km s“^ with a 
standard deviation of 7.6 km s“^ for the evolved stars. Note 
that our results are given for solar abundances, but even 
though we find slightly lower values of microturbulence ve¬ 
locity (by nearly 2 km s“^) when we use the abundances 
derived by the various authors, our values are still higher 
than those from the literature. 

Fixing or constraining the microturbulence velocity 
with only a few lines can lead to large errors for the 
other parameter estimates as well. Figure shows an 
example of the most probable pairs [Te//, log(<?)l for 
5 = 2, 6 , 10, 14, and 18 km s“^ when calculating the 
likelihoods with the lines He I A4388, He I±Mg II around 
4476 A, Sill AA4128,4130, and Si III AA4552,4567,4574 


of the artificial spectrum with the following specifica¬ 
tions: Te// = 21350 K, log((;) = 3.8, 5 = 3 km s“'^, and 
V sin(i) = 15 km s“^. This shows that each line has a 
different sensitivity to the microturbulence velocity and, 
depending on the line or the group of lines chosen as a 
diagnostic, the result differs in Te// and log((;) but also in 
5 . For instance here, the most probable value of 5 given by 
the helium lines is 10 km s“^ (when marginalizing over all 
the other parameters) while the silicon lines give 2 km s“^. 
Thus, it is very important to consider 5 as a free parameter 
as is usually done with Te//, log((;), and v sin(i), and to 
constrain these four parameters simultaneously with as 
many lines as possible. Also it is even more important to 
consider the microturbulence velocity as a free parameter 
when studying the abundance, since it has a significant 
impact on this parameter as well (Nieva fc Przybilla|2010 1 . 


5 CONCLUSIONS 

We described here a spectral analysis method based on 
Bayesian statistics that simultaneously constrains four stel¬ 
lar parameters (effective temperature, surface gravity, pro- 
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Figure 21. The most probable pair log(g)] as a function 

of ^ for different lines (as indicated in the plot) for the same arti¬ 
ficial spectrum as in Figure The behaviour as a function of g 
is shown in the plane log(g)]. The likelihoods are marginal¬ 

ized over all v sin(i) values. 


jected rotational velocity, and microturbulence velocity). 
The method cumulates the information provided simultane¬ 
ously by all available lines in a given spectrum. This allow 
us to find the best global solution for all the lines and all 
the parameters at the same time, while providing reduced 
uncertainties. The Bayesian formalism naturally gives the 
uncertainties for each parameter in relation with the other 
parameter possible values, with the data uncertainties, and 
with the limitations of the models used. 

This method is completely objective since it does not 
rely on the user’s visual judgement, nor does it depend on 
starting-parameter estimates. We have shown that adopting 
an inaccurate starting value or fixing the value of a given 
parameter, along with using a classical iterative method, 
can lead to substantial errors on the other parameters. Also, 
because our method considers all the parameters and all the 
lines at the same time, it does not depend on the sensitivity 
of a specific diagnostic line. 

This method is also self-consistent and efficient even 
when it is applied to spectra with high noise level or when 
underestimating the quality of the data. Moreover, since our 
method cumulates the constraints of all available lines, it is 
still able to give accurate results with small uncertainties 
even though nearly all the lines are heavily blended due to 
the modest resolution (2.3 A) of our data. 

Note also, that the method actually works with atmo¬ 
spheric models from TLUSTY but also from ATLAS9 and 
PHOENIX dBaron, Chen fc HauschiIdt|[Mol ), and can be 
easily adapted to virtually any atmospheric model. Further¬ 
more, it is a fast method, once the basic preparation is done 
(preparation of the data as suitable input for the code and 
creation of the basic grids) the complete basic parameters 
analysis of a star takes less than 3 minutes on a typical 
portable computer (we use a Pentium(R) Dual-Core CPU 
T4300 @ 2.10GHz with 4 Go of RAM). 

Finally, the comparison of our results with those from 


the literature is very satisfying overall. Most of the differ¬ 
ences found are easily explained in terms of the method used 
(an iterative method versus our simultaneous method) and 
the number and choice of the diagnostic lines used (few lines 
for specific parameters versus, as we did, all available lines 
with the same weight for all the parameters). One important 
behaviour is pointed out in this study: the microturbulence 
velocity is often underestimated for the B dwarfs as well as 
for the giants and supergiants. Furthermore, we also con¬ 
firmed that cluster B stars are on average faster rotators 
than field B stars. 

Now that we have demonstrated that we can success¬ 
fully apply our method to gather stellar parameters using 
solar abundances, we will investigate in a future paper the 
impact and efficiency of our method for the abundance de¬ 
termination of various chemical elements. 
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APPENDIX A: SUPPLEMENTARY FIGURES 



Figure Al. Illustration of the commutativity of the analysis. The lines from Fig are used here with a similar grid of synthetic 
spectra where all the parameters other than T^ff and ^ have been marginalized. The green dash-dotted frame on the left shows each 
line considered with its respective likelihood, while the blue dashed frames holds the posterior probability distribution obtained when 
considering several possible combinations of the likelihoods. From the top to the bottom, each column of the blue frame shows the 
evolution of the posterior probability distribution as new lines are added to the analysis. While the order in which the lines are added 
differs from column to column, the final posterior probability distribution (bottom row) remains the same in each case. Note that the 
products written in each graphic (such as L1*L2*L3) are simplified for clarity. In fact, each posterior distribution is the normalised 
product of the likelihoods as stated in equation]^ 
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APPENDIX B: TABLES OF STELLAR PARAMETERS 


Table Bl: Basic Parameters of Field B Stars 


Star Id. 

Teff 

log(5) 

V sin(i) 

i 

ref 


(K) 


(km s“^) 

(km s 


HD 886 

21250 ± 130 

3.98 ± 0.04 

31 

8 ± 1 

TW 


21600 ± 1000 

3.7 ± 0.3 

8± 2 

3.3 

AND99 


23000 ± 1000 

3.8 ± 0.1 

10 ± 1 

1 

LEFIO 


21667 ± 650 

3.8 ± 0.2 

0± 9 

3± 2 

TAKIO 

HD3827 

29440 talo 

4.27 ± 0.04 

155 ± 5 

6 ± 1 

TW 


24678 ±450 

3.836 ± 0.038 

164 ± 10 


HGIO 

HD5882 

18250 ± 130 

4.27 ± 0.04 

330 

7 +0 
' -2 

TW 


15626 ± 150 

3.731 ± 0.021 

358 ± 5 

2 

HGIO 

HD32249 

20000 ± 130 

4.26 ± 0.13 

90 ± 10 

5 ± 1 

TW 


17081 ± 150 

3.859 ± 0.025 

60 ± 12 

2 

HGIO 


18300 ± 300 

3.9 ± 0.1 

46 

0 6 
+ 1 

LYU05 


18500 ± 1000 

4. ± 0.3 

48 ± 8 

4 

AND99 


18890 ± 567 

4.1 ± 0.2 

30 ±9 

3± 2 

TAKIO 

HD32672 

20950 tiso 

3.93 ± 0.04 

37 ± 13 

4 la 

TW 


18089 ± 300 

3.434 ± 0.034 

32 ± 14 

2 

HGIO 

HD33816 

22430 tffo 

3.94 ± 0.06 

0 +12 
° -8 

7 ± 1 

TW 


21613 ± 500 

3.818 ± 0.049 

24 ± 24 


HGIO 

HD35912 

19400 tigg 

A QQ +0.08 
_0.04 

1 5 

20 -15 

1 1 ? 

TW 


19000 ± 1000 

4.1 ± 0.3 

22 ± 4 

4 

AND99 

HD36280 

28000 ± 130 

4.02 ± 0.04 

203 

9 tl 

TW 


23611 ± 600 

3.613 ± 0.053 

206 ± 9 

2 

HGIO 

HD37032 

30000 t° 5 o 

4.07 ± 0.04 

63 ± 15 

10 1 ^ 

TW 


22499 ± 300 

3.442 ± 0.028 

70 ± 14 

2 

HGIO 

HD40160 

20010 ^240 

4 05 +0 07 
+-UO -0.05 

128 tl 

5 I 2 

TW 


16453 ± 250 

3.735 ± 0.027 

115 ± 13 


HGIO 

HD42782 

15000 

4.02 ± 0.04 

60 til 

0 

O’ -0 

TW 


15133 ± 150 

3.790 ± 0.023 

44 ± 21 

2 

HGIO 

HD44700 

17380 ± 130 

4.24 ± 0.04 

98 +1* 

^9 -22 

5 ± 1 

TW 


17520 

4.24 ± 0.13 

/ +10 

4 -4 

6± 1 

TW^ 


15190 ± 150 

3.907 ± 0.026 

26 ± 26 

2 

HGIO 


17500 ± 500 

3.8 ± 0.1 

8± 1 

6± 2 

LEFIO 


16551 ±497 

4.2 ± 0.2 

0± 9 

3± 2 

TAKIO 

HD45418 

16750 ± 130 

4.27 ± 0.04 

247 t? 

6 ± 1 

TW 


16000 ± 500 

4. ± 0.2 

237 ± 28 

3± 2 

LEFIO 

HD77770 

19250 ± 130 

3.77 ± 0.04 

rn _j4 

7 ± 1 

TW 


19000 ± 130 

3.53 ± 0.04 

O' -0 

4 I 2 

TW^ 


16935 ± 200 

3.386 ± 0.023 

49 ± 49 

2 

HGIO 

HD160762 

17050 ± 200 

3.83 ± 0.03 

0 +12 
■0 -3 

8 ± 1 

TW 


17000 ± 250 

3.77 ± 0.12 

15 

10+1 
-'•9 -1.3 

LYU05 


19500 ± 500 

4.1 ± 0.1 

8± 2 

3± 2 

LEFIO 


17440 ± 523 

3.9 ± 0.2 

0± 9 

3± 2 

TAKIO 

HD166182 

19130 

3.54 ± 0.04 

0 +00 

O’ -0 

10 ± 1 

TW 


20000 ± 600 

3.54 ± 0.16 

44 

4-7 ti. 

LYU05 

HD184171 

16980 

0 on +0.05 
-0.07 

0 +7 

■2 -3 

3 ± 1 

TW 


16100 ± 700 

3.62 ± 0.15 

29 

0 ± 0 

LYU05 


15858 ±476 

3.5 ± 0.2 

15 ± 9 

3± 2 

TAKIO 

Continued on next page 
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Table B1 - 

continued from 

previous page 


star Id. 

Teff 

log(5) 

V sin(i) 

c 

ref 


(K) 


(km s“^) 

(km s 


HD185330 

15000 “ 

3.70 ± 0.04 

115 tm 

1 ± 1 

TW 


16167 ± 485 

3.7 ± 0.2 

15 ± 9 

3± 2 

TAKIO 

HD188892 

15120 

3.57 ± 0.04 

-3 

4 ± 1 

TW 


I 4 OO 8 ± 420 

3.3 ± 0.2 

30 ± 9 

3± 2 

TAKIO 

HD189944 

15500 ± 130 

4.02 ± 0.04 

-0 

6 ± 1 

TW 


14793 ± 444 

4.0 ± 0.2 

10 ± 9 

3± 2 

TAKIO 

HD190066 

23560 tiso 

2.90 ± 0.04 

99 

25 

TW 


21000 ± 1000 

2.88 ± 0.25 

82 ± 20 

15 

SEA08 

HD191243 

15130 

2.73 ± 0.04 

-3 

11 ± 1 

TW 


13923 ± 4 I 8 

2.5 ± 0.2 

30 ± 9 

3± 2 

TAKIO 


14500 ± 1000 

2.75 ± 0.25 

38 ± 20 

20 

SEA08 


I 48 OO ± 1000 

2.6 ± 0.2 

38 ± 10 

8± 5 

MAR08 

HD192660 

25250 ± 130 

2.90 ± 0.04 

145 ± 5 

28 ± 1 

TW 


30000 ± 1000 

3.25 ± 0.25 

94 ± 20 

20 

SEA08 

HD193183 

20480 

cy yy +0.04 
* ‘ -0.08 

44 « 

21 ± 1 

TW 


18500 ± 1000 

2.63 ± 0.25 

68 ± 20 

20 

SEA08 

HD196035 

18290 

4.20 ± 0.06 

0 +12 

■2 -3 

8 ± 1 

TW 


17100 ± 250 

4.0 ± 0.08 

34 

2± 2 

LYU05 


17499 ± 525 

4.3 ± 0.2 

20 ± 9 

3± 2 

TAKIO 

HD198478 

16380 tiao 

2.09 ± 0.04 

u -0 

31 ± 1 

TW 


17500 ± 1000 

2.25 ± 0.25 

61 ± 20 

20 

SEA08 


17500 ± 1000 

2.1 ± 0.2 

39 ± 10 

8± 5 

MAR08 

HD198820 

16760 tiao 

3.77 ± 0.04 

-3 

6 ± 1 

TW 


15600 ± 300 

3.65 ± 0.15 

35 

n 1 

- 0.1 

LYU05 


15852 ±476 

3.8 ± 0.2 

15 ± 9 

3± 2 

TAKIO 

HD199578 

15000 

4.05 ± 0.04 

-3 

5 ± 1 

TW 


I 448 O ± 434 

4.0 ± 0.2 

30 ± 9 

3± 2 

TAKIO 

HD202347 

22630 tiao 

4.08 ± 0.04 

95 ± 5 

8 ± 1 

TW 


23280 ± 900 

4.13 ± 0.1 

119 ± 6 

6 ± 1.5 

DAF07 

HD205794 

27980 

4.17 ± 0.03 

c +25 

0 -5 

9 ± 1 

TW 


26890 ± 1000 

4.21 ± 0.1 

14 ± 3 

8 ± 1.5 

DAF07 

HD206165 

19020 ti 5 o 

2.59 ± 0.04 

75 ± 5 

21 ± 1 

TW 


18000 ± 1000 

2.50 ± 0.25 

73 ± 20 

15 

SEA08 


19300 ± 1000 

2.5 ± 0.2 

45 ± 10 

8± 5 

MAR08 

HD206540 

15000 “ 

4.21 ± 0.04 

0 -6 

5 ± 1 

TW 


13981 ± 419 

4.0 ± 0.2 

20 ± 9 

3± 2 

TAKIO 

HD209419 

15020 tfo° 

3 96 +0 0 ^ 

_o.o5 

40 t\t 

5 ± 1 

TW 


14100 ± 600 

3.70 ± 0.12 

17 

0± 0 

LYU05 


14404 ± 432 

3.8 ± 0.2 

20 ± 9 

3± 2 

TAKIO 

HD213087 

24870 tiao 

2.90 ± 0.04 

103 ± 8 

24 ± 1 

TW 


27000 ± 1000 

3.13 ± 0.25 

88 ± 20 

20 

SEA08 

HD214263 

20860 ^240 

4.15 ± 0.04 

90 


TW 


19027± 250 

3.807 ± 0.025 

63 ± 16 

2 

HGIO 

HD214432 

19500 ± 130 

-0.05 

135 ± 5 

7 ± 1 

TW 


16124 ± 200 

3.890 ± 0.028 

137 ± 9 

2 

HGIO 

HD215191 

22010 tlfo 

4.0 ± 0.06 

200 ± 5 

8 ± 1 

TW 


21487 ± 350 

3.83 ± 0.03 

195 ± 11 

2 

HGIO 

Continued on next page 
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Table B1 - continued from previous page 


Star Id. 

Teff 

(K) 

log(s) 

V sin(f) 

(km s“^) 

C 

(km s 

ref 

HD215371 

20380 tiaE; 
17688 ± 250 

4.27 ± 0.04 
3.806 ± 0.024 

165 ± 5 

138 ± 10 

8 ± 1 

2 

TW 

HGIO 



1: Spectrum with an effective resolution of 0.9 A. 

Referenc es: TW = This Work using mod els from T LUSTY with the sola r composit ion from Asplund et al. (20091; 

(|2010|; TAKIO = |Takeda et al.| dlSOlOl); 


HGIO = 
SEA08 = 


Huang, Gies fc McSwain (20101; LEFIO = 


HG06 = Huang & Gies 


Searie et al. (200^; MAR08 = 


Lefever et al. 


Markova & Puls 


(2008 




DAF07 = 


'{ 2006a b I; LYU05 = Lyubimkov et al. (20051; and AjND 99 = Andrievsky et al. (19991. 


Daflon 


et al. 


I 20071 
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Table B2: Basic Parameters of Cluster B Stars 


Cluster Id. 

Star Id. 

Teff 

log(5) 

V sin(i) 

C 

ref 



(K) 


(km s“^) 

(km s ") 


NGC1960 

BD+34 1103s 

19660 

4.30 ± 0.04 

288 12 

4 11 

TW 



11291 ± 150 

3.983 ± 0.015 

11 ±1 

a 

HGIO 


SAB20 

15000 

4 20 

-0.08 

301 1“ 

6 ± 3 

TW 



14845 ± 100 

4.242 ± 0.018 

141 ± 1 

2 

HGIO 


SAB22 

15000 4o“ 

4.35 ± 0.04 

118 12 

3 11 

TW 



13519 ± 100 

4.011 ± 0.019 

135 ± 6 

2 

HGIO 


SAB23 

15000 

-0.06 

300 ± 20 

^ -0 

TW 



13115 ± 150 

4.261 ± 0.028 

302 ± 11 


HGIO 


SAB30 

15240 +_IZ 

4 44 +0-16 

351 1« 

IQ +8 

-13 

TW 



13651 ± 150 

4.01 ± 0.03 

19 ± 24 

2 

HGIO 


SA058227 

16870 

4.47 ± 0.03 

87 11" 

2 11 

TW 



14625 ± 100 

4.031 ± 0.011 

12 ±1 

2 

HGIO 


SAO58230 

20830 

4.09 ± 0.04 

95 

_;LQ 

6 ± 1 

TW 



19611 ± 300 

3.819 ± 0.026 

57 ± 11 

a 

HGIO 


SA058239 

20600 

4.06 ± 0.04 

315 ± 5 

5 ± 1 

TW 



19582 ± 200 

3.911 ± 0.018 

199 ± 1 

2 

HGIO 

NGC884 

BD+56 556 

25130 tlt°o 

0 qy +0.04 
-0.08 

140 ± 5 

10 ± 1 

TW 



25931 ±401 

4.160 ± 0.041 

69 ±9 

a 

HG06 


BD+56 604 

23250 ± 130 

3.83 ± 0.04 

136 12 

4 ± 1 

TW 



23912 ± 500 

3.999 ± 0.048 

91 ± 10 

2 

HGIO 


HD14321 

22620 1122 

3.77 ± 0.04 

150 ± 5 

6 ± 1 

TW 



23601 ± 550 

3.896 ± 0.046 

106 ± 10 

2 

HGIO 


HD 14476 

24630 1122 

3.45 ± 0.04 

105 ± 5 

13 ± 1 

TW 



26006 ± 800 

3.149 ± 0.094 

40 ± 14 

2 

HGIO 


HG2048 

19720 1112 

4.31 ± 0.04 

203 12 

6 11 

TW 



18001 ± 198 

3.944 ± 0.021 

125 ± 8 

2 

HG06 


HG2255 

21670 1122 

3.83 ± 0.04 

329 12 

6 ± 1 

TW 



20360 ± 118 

3.560 ± 0.022 

318 ± 10 

2 

HG06 


HG2794 

22510 ± 260 

q 70 +0.07 

0 . (0 _o .05 

140 ± 10 

8 ± 1 

TW 



21624 ± 291 

3.611 ± 0.021 

19 ± 8 

2 

HG06 


W168 

21370 1122 

4.27 ± 0.04 

29 112 

5 

^ -2 

TW 



19259 ± 350 

3.98 ± 0.04 

44 ± 14 

2 

HGIO 


HGIO = 


Huang, Gies & McSwain 

(2010 

1 ; and HG06 = 

Huang & Gies 

(2006a b 
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Table B3: Basic Parameters of Field B Stars with Modified He and Si 
Abundances from Lefever et al. (20101 


Star Id. 

Teff 

(K) 

log(ff) 

V sin(z) 
(km s~^) 

C 

(km s 

log(^) + 12 

log(f) + 12 

ref 

HD160762 

18090 

3.97 ± 0.04 

O +7 
-3 

9 ± 1 

11 

7.14 

TW 


19500 ± 500 

4.1 ± 0.1 

8 ± 2 

3± 2 



LEFIO 

HD44700 

17250 ± 130 

4.22 ± 0.04 

S +16 
® -8 

6 ± 1 

11 

7.07 

TW 


17500 ± 500 

3.8 ± 0.1 

8 ± 1 

6 ± 2 



LEFIO 

HD45418 

16750 ± 130 

4.27 ± 0.04 

247 

6 ± 1 

11 

7.51 

TW 


16000 ± 500 

4 . ± 0.2 

237 ± 28 

3± 2 



LEFIO 

HD 886 

21120 tlfo 

3.95 ± 0.04 

34 ± 8 

9 ± 1 

11 

7.16 

TW 


23000 ± 1000 

3.8 ± 0.1 

10 ± 1 

1 



LEFIO 


References: TW = This Work; LEFIO = 


Lefever et al. 


( 20101 . 
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Table B4: Basic Parameters of Field B Stars with Modified O and Ne 
Abundances from Takeda et al. (20101 


Star Id. 

Teff 

log(3) 

V sin(i) 


log(g) + 12 

log(^) + 12 

ref 


(K) 


(km s“^) 

(km s 




HD160762 

16880 ± 130 

3.81 ± 0.04 

q +12 

^ -3 

8l2 

8.74 

8.08 

TW 


17440 ± 523 

3.9 ± 0.2 

o± 9 

3± 2 



TAKIO 

HD184171 

17000 ± 130 

3.86 ± 0.04 

0 +7 
-3 

3 ± 1 

8.69 

8.01 

TW 


15858 ±476 

3.5 ± 0.2 

15 ± 9 

3± 2 



TAKIO 

HD185330 

15000 

3.70 ± 0.04 

120 ± 12 

2 tl 

8.35 

8.15 

TW 


16167± 485 

5.7± 0.2 

15 ± 9 

3± 2 



TAKIO 

HD188892 

15120 t\fo 

3.57 ± 0.04 

0 +7 
-3 

4 ± 1 

8.67 

8.01 

TW 


I4OO8 ± 420 

3.3 ± 0.2 

30 ± 9 

3± 2 



TAKIO 

HD189944 

15500 ± 130 

4.02 ± 0.04 

0 

^ -0 

6 ± 1 

8.79 

8.11 

TW 


14793 ± 444 

4.0 ± 0.2 

10± 9 

3± 2 



TAKIO 

HD191243 

15230 till 

9 '7^ +0.05 
^+0 -0.07 

0 +7 
^ -3 

10 

8.66 

7.95 

TW 


13923 ± 4I8 

2.5 ± 0.2 

30 ± 9 

3± 2 



TAKIO 

HD196035 

18350 

4.20 ± 0.03 

q +12 

^ -3 

8-1 

8.7 

8.08 

TW 


17499 ± 525 

4.3 ± 0.2 

20 ± 9 

3± 2 



TAKIO 

HD198820 

16750 

3.77 ± 0.04 

0 +7 
^ -3 

6 ± 1 

8.8 

8.13 

TW 


15852 ±476 

3.8 ± 0.2 

15 ± 9 

3± 2 



TAKIO 

HD199578 

15000 

4.05 ± 0.04 

0 +7 
^ -3 

Gtl 

8.7 

8.14 

TW 


I448O ± 434 

4.0 ± 0.2 

30 ± 9 

3± 2 



TAKIO 

HD206540 

15000 tf° 

4 21 

-0.05 

4 2 ® 

4 ± 1 

8.72 

8 

TW 


13981 ± 419 

4.0 ± 0.2 

20 ± 9 

3± 2 



TAKIO 

HD209419 

15000 

3.95 ± 0.04 

41 11 ; 

5 ± 1 

8.73 

7.99 

TW 


14404 ± 432 

3.8 ± 0.2 

20 ± 9 

3± 2 



TAKIO 

HD32249 

19980 

4.26 ± 0.13 

86 ± 8 

5 ± 1 

8.69 

7.95 

TW 


18890 ± 567 

4.1 ± 0.2 

30 ± 9 

3± 2 



TAKIO 

HD44700 

17370 

4.24 ± 0.04 

29 ± 10 

5 ± 1 

8.73 

8.04 

TW 


16551 ±497 

4.2 ± 0.2 

0± 9 

3± 2 



TAKIO 

HD886 

20930 

3.95 ± 0.04 

31 

8 ± 1 

8.77 

8.09 

TW 


21667 ± 650 

3.8 ± 0.2 

0± 9 

3± 2 



TAKIO 


References: TW = This Work; TAKIO = |Takeda et al. ( 2010| ). 
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Table B5: Basic Parameters of Field B Stars with Modified C, N, O 
Abundances from Searle et al. (20081 


Star Id. 

Teff 

(K) 

log(5) 

V sin(i) 

(km s“^) 

e 

(km s 

log(S) + 12 

log(f) + 12 

log(§) + 12 

ref 

HD 190066 

23750 ± 190 

2.90 ± 0.04 

95 ± 9 

25 ± 1 

7.88 

8.15 

8.53 

TW 


21000 ± 1000 

2.88 ± 0.25 

82 ± 20 

15 




SEA08 

HD191243 

15120 

2.70 ± 0.04 

o +7 
-3 

11 ± 1 

7.70 

7.65 


TW 


14500 ± 1000 

2.75 ± 0.25 

38 ± 20 

20 




SEA08 

HD192660 

25140 ± 160 

2.88 ± 0.03 

140 ± 5 

29 71 

8.02 

7.51 

8.73 

TW 


30000 ± 1000 

3.25 ± 0.25 

94 ± 20 

20 




SEA08 

HD193183 

19610 +_l% 

2.65 ± 0.04 

32 ± 12 

23 ± 1 

7.66 

8.15 

8.73 

TW 


18500 ± 1000 

2.63 ± 0.25 

68 ± 20 

20 




SEA08 

HD 198478 

16130 tlfo 

2.04 ± 0.04 

^ -0 

33 ± 1 

7.86 

8.29 

8.45 

TW 


17500 ± 1000 

2.25 ± 0.25 

61 ± 20 

20 




SEA08 

HD206165 

18630 71^ 

2.55 ± 0.04 

75 ± 5 

23 ± 1 

7.96 

8.15 

8.43 

TW 


18000 ± 1000 

2.50 ± 0.25 

73 ± 20 

15 




SEA08 

HD213087 

24620 7112 

2.86 ± 0.04 

103 ± 8 

24 ± 1 

8.00 

8.15 

8.73 

TW 


27000 ± 1000 

3.13 ± 0.25 

88 ± 20 

20 




SEA08 


References; TW = This Work; SEA08 = 


Searle et al. 


(2008||. 


Table B6: Basic Parameters of Field B Stars with Modified He and Si 
Abundances from [Markova fc Puls| (|2008|) 


Star Id. 

Tff 

(K) 

log(5) 

V sin(i) 
(km s“^) 

7 

(km s 

log(l^) + 12 

log(f) + 12 

ref 

HD191243 

15240 7112 

Q 'T'T +0.04 
^ * -0.08 

O +7 
-3 

11 ± 1 

10.954 

7.48 

TW 


14800 ± 1000 

2.6 ± 0.2 

38 ± 10 

8± 5 



MAR08 

HD 198478 

16270 7112 

to 

b 

00 

1 + 
o o 
b b 
00 

o +28 
^ -2 

26 71 

11.176 

7.58 

TW 


17500 ± 1000 

2.1 ± 0.2 

39 ± 10 

8± 5 



MAR08 

HD206165 

19380 7112 

2.63 ± 0.04 

76 72 

16 7? 

11.176 

7.58 

TW 


19300 ± 1000 

2.5 ± 0.2 

45 ± 10 

8± 5 



MAR08 


(|2008l. 


References: TW = This Work; MAR08 = 


Markova & Puls 
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Table B7: Basic Parameters of Field B Stars with Modified C, N, O, Mg, 
Al, and Si Abundances from Daflon et al. (20071 


Star Id. 

Rff 

(K) 

log ( 5 ) 

V sin(i) 
(km s~^) 

(km s 

<C) e{N) 

e(0) 

e{Mg) 

RAl) 

e{Si) 

ref 

hd202347 

23730 tlto 
23280 ± 900 

4 27 

-0.08 

11 ± 0.1 

96 

119 ± 6 

9 ± 1 

6 

0 -1.5 

7.54 

8.54 

7.35 

6.28 

7.13 

TW 

DAF07 

hd205794 

25010 ligj; 
26890 ± 1000 

3.95 ± 0.04 

4.2 ± 0.1 

q +1T 

0 -3 

14 ± 3 

11 ± 1 

° -1.5 

8.09 7.57 

8.63 



7.09 

TW 

DAF07 


References: TW = This Work; AND99 = ^Daflon et al. (20071. Here, e(A’) 


log(X/H) + 12. 


Table B8: Basic Parameters of Field B Stars with Modified He and Mg 
Abundances from [Lyubimkov et ah] (|2005| 


Star Id. 

Teff 

(K) 

log ( 5 ) 

V sin(i) 

(km s-i) 

(km s-i) 

log(ll^) + 12 

log(^) + 12 

ref 

HD160762 

17280 

3.88 ± 0.04 


7 ± 1 

11.057 

7.62 

TW 


17000 ± 250 

5.77± 0.12 

15 

y 0 +1 
-1.3 



LYU05 

HD166182 

19150 

0 rc +0.07 

0.00 _o.05 

n + 1 ° 

-0 

9 ± 1 

11.045 

7.39 

TW 


20000 ± 600 

3.54 ± 0.16 

44 

4.7 tl. 



LYU05 

HD184171 

16880 

3.85 ± 0.04 

0 +7 
-3 

2 ± 1 

11.041 

7.89 

TW 


16100 ± 700 

3.62 ± 0.15 

29 

o± 0 



LYU05 

HD196035 

18210 till 

4.18 ± 0.04 

0 +7 
-3 

6 

11.065 

7.66 

TW 


17100 ± 250 

4.0 ± 0.08 

34 

2± 2 



LYU05 

HD198820 

16500 ± 130 

3.73 ± 0.04 


3 ± 1 

11.11 

8.01 

TW 


15600 ± 300 

3.65 ± 0.15 

35 

0 1 

-0.1 



LYU05 

HD209419 

15000 

3.95 ± 0.04 

60 ± 5 

4 ± 1 

11.093 

8.03 

TW 


14100 ± 600 

3.70 ± 0.12 

17 

o± 0 



LYU05 

HD32249 

18880 

4.11 ± 0.04 

90 ± 5 

2 ± 1 

11.124 

8.02 

TW 


18300 ± 300 

3.9 ± 0.1 

46 

/7 +0 

-0.5 



LYU05 


References: TW = This Work; LYU05 = 


Lyubimkov et al. (20051. 
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Table B9: Basic Parameters of Field B Stars with Modified C and N 
Abundances from Andrievsky et al. (19991 


Star Id. 

Tff 

(K) 

log(fl) 

V sin(i) 
(km s“^) 

(km s 

log(S) + 12 

log(f) + 12 

ref 

HD32249 

20000 ± 130 

4.26 ± 0.13 

86 

5 ± 1 

8.19 

7.8 

TW 


18500 ± 1000 

4. ± 0.3 

48 ± 8 

4 



AND 99 

HD35912 

19610 

4.27 ± 0.04 

17 

J- i _i7 

1 

8.28 

7.76 

TW 


19000 ± 1000 

4.1 ± 0.3 

22 ± 4 




AND 99 

HD886 

21250 ± 130 

3.97 ± 0.04 

4 If 

8 ± 1 

8.11 

7.55 

TW 


21600 ± 1000 

3.7 ± 0.3 

8± 2 

3.3 



AND 9 9 


References: TW = This Work; AND99 = [Andrievsky et al. (19991. 


Table BIO: Basic Parameters of Field B Stars Returned by ATLAS9 with 
Solar Composition from|Grevesse fc Sauval| ||1998[) 


Star Id. 

Teff 

log(ff) 

V sin(i) 


ref 


(K) 


(km s“^) 

(km s 


HD214263 

19760 Ifo 

4.08 ± 0.04 

88 

It? 

TW 

(16) 

18570 ± 350 

3.87 ± 0.06 

63 

2 

TWHG 


19027± 250 

3.807 ± 0.025 

63 ± 16 

2 

HGIO 

HD214432 

18870 Ifo 

4.47 ± 0.04 

132 

5 ± 1 

TW 

(17) 

15560 

3.95 ± 0.08 

137 

2 

TWHG 


16124 ± 200 

3.890 ± 0.028 

137 ± 9 

2 

HGIO 

HD215191 

20500 

3.91 ± 0.04 

200 ± 5 

6 ± 1 

TW 

(18) 

20490 ± 360 

3.88 ± 0.06 

195 

2 

TWHG 


21487 ± 350 

3.83 ± 0.03 

195 ± 11 

2 

HGIO 

HD215371 

19380 

4.23 ± 0.04 

160 ± 5 

5 ± 1 

TW 

(19) 

18131 ± 450 

3.99 ± 0.1 

138 

2 

TWHG 


17688 ± 250 

3.806 ± 0.024 

138 ± 10 

2 

HGIO 

HD32249 

18890 ± 130 

4.25 ± 0.04 

95 ± 5 

4 ± 1 

TW 

(20) 

16490 ± 450 

3.84 ± 0.11 

60 

2 

TWHG 


17081 ± 150 

3.859 ± 0.025 

60 ± 12 

2 

HGIO 

HD32672 

19770 ttTo 

0 07 +0.07 
-0.05 

18 Ifs 

1 ± 1 

TW 

(21) 

18290 ± 570 

3.69 ± 0.08 

32 

2 

TWHG 


18089 ± 300 

3.434 ± 0.034 

32 ± 14 

2 

HGIO 

HD33816 

20780 Ilf 

3.85 ± 0.04 

8 If 

4 ± 1 

TW 

(22) 

21110 ± 1540 

3.91 ± 0.06 

24 

2 

TWHG 


21613 ± 500 

3.818 ± 0.049 

24 ± 24 

2 

HGIO 

HD36280 

24330 Ilf 

3.91 ± 0.04 

214 ± 6 

8 ± 1 

TW 

(23) 

24020 ± 530 

3.84± 0.1 

206 

2 

TWHG 


23611 ± 600 

3.613 ± 0.053 

206 ± 9 

2 

HGIO 

HD37032 

29540 tlto 

A 98 +0.07 

4.Z5 _o,05 

43 +ll 

5 ± 1 

TW 

(24) 

22110 ± 1300 

3.56 ± 0.17 

70 

2 

TWHG 


22499 ± 300 

3.442 ± 0.028 

70 ± 14 

2 

HGIO 

HD3827 

29500 ± 500 

4.43 ± 0.04 

165 ± 5 

6 ± 1 

TW 

(25) 

24100 ± 810 

3.93 ± 0.14 

164 

2 

TWHG 


24 678 ± 450 

3.836 ± 0.038 

164 ± 10 

2 

HGIO 

Continued on next page 
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Table BIO 

- continued from 

previous page 



Star Id. 

Teff 

log(ff) 

V sin(i) 


ref 


(K) 


(km s“^) (k 

ms 


HD42782 

13040 

0 79 +0.07 

0 ^ (^ -0.05 

155 

0 

1 + 

0 to 

TW 

(26) 

14570 ± 590 

3.9 ± 0.16 

44 

2 

TWHG 


15133 ± 150 

3.790 ± 0.023 

44 ± 21 

2 

HGIO 

HD44700 

16900 

00 

0 0 

d d 

+ 1 

CO 

CM 

82 

3 

TW 

(27) 

15010 ± 510 

3.91 ± 0.12 

26 

2 

TWHG 


15190 ± 150 

3.907 ± 0.026 

26 ± 26 

2 

HGIO 

HD5882 

17500 ± 130 

4.27 ± 0.04 

328 ± 6 

6 ± 3 

TW 

(28) 

14760 ± 300 

3.71 ± 0.09 

358 

2 

TWHG 


15626 ± 150 

3.731 ± 0.021 

358 ± 5 

2 

HGIO 

HD77770 

18500 ± 130 

3.78 ± 0.04 

1 c +5 
-10 

1 -H 0 
+ 1 

0 

TW 

(29) 

15590 ±420 

3.31 ± 0.07 

49 

.8 

TWHG 


16935 ± 200 

3.386 ± 0.023 

49 ± 49 

2 

HGIO 


References: TW = This Work; T WHG = This Work using H 7 only; 
HGIO = 


Huang, Gies & McSwain (20101. 

Numbers between parenthesis refer to the ’Star id n°’ in Figure 
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Table Bll: Stellar Parameters of Gluster B Stars Returned by ATLAS9 
with Solar Composition from Grevesse & Sauval (1998 1 


Cluster Id. 

Star Id. 

Rff 

log( 5 ) 

V sin(i) 

i 

ref 



(K) 


(km s~^) 

(km s 


NGC1960 

BD341103S 

18870 

4.31 ± 0.04 

290 Ig 

2 ± 1 

TW 


( 0 ) 

20030 ± 920 

4.40 ± 0.17 

77 

2 

TWHG 



17297± 150 

3.983 ± 0.015 

77 ±7 

2 

HGIO 


SAB20 

12970 

Q O'? +0.05 
'J-o' -0.07 

274 IS 

1 

TW 


( 8 ) 

14560 ± 620 

4.25 ± 0.15 

141 

2 

TWHG 



14845 ± 100 

4.242 ± 0.018 

141 ± 7 

2 

HGIO 


SAB22 

14500 ± 130 

4.33 ± 0.04 

126 IS 

3 ± 1 

TW 


(9) 

13500 ±210 

4.09 ± 0.06 

135 

2 

TWHG 



13579 ± 100 

4.071 ± 0.019 

135 ± 6 

2 

HGIO 


SAB23 

12780 

4 07 

-0.08 

244 IS 

4 ± 2 

TW 


( 10 ) 

12910 ± 1220 

4.18 ± 0.3 

302 

2 

TWHG 



13715 ± 150 

4.267 ± 0.028 

302 ± 11 

2 

HGIO 


SAB30 

11920 IS!; 

4.89 

261 +_fo 

12 ± 3 

TW 


( 11 ) 

12160 ± 840 

4.88 ± 0.3 

79 

2 

TWHG 



13651 ± 150 

4.07 ± 0.03 

79 ± 24 

2 

HGIO 


SA058227 

16600 tlfo 

-0.07 

123 IS 

5 ± 1 

TW 


( 12 ) 

14790 ± 250 

4.19 ± 0.08 

72 

2 

TWHG 



14625 ± 100 

4.031 ± 0.017 

72 ±7 

2 

HGIO 


SAO58230 

20190 

4.11 ± 0.04 

98 ± 12 

2 

TW 


(13) 

19420 ± 630 

4.05 ± 0.1 

57 

2 

TWHG 



19671 ± 300 

3.879 ± 0.026 

57 ± 11 

2 

HGIO 


SA058239 

19120 tlfo 

3.95 ± 0.04 

315 ± 5 

2 ± 1 

TW 


(14) 

21000 ± 560 

4.15 ± 0.08 

199 

2 

TWHG 



19582 ± 200 

3.911 ± 0.018 

199 ± 7 

2 

HGIO 

Continued on next page 
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Table Bll - continued from previous page 


Cluster Id. 

Star Id. 

Teff 

log(5) 

V sin(i) 


ref 



(K) 


(km s“^) 

(km s 

NGG884 

BD56556 

24390 1 1 ^ 5 ; 

0 Q'T +0.07 
-0.05 

145 ± 5 

9 ± 1 

TW 


( 1 ) 

25000 ± 2300 

4.25 ± 0.25 

69 

2 

TWHG 



25937 ± 401 

4.160± 0.041 

69 ± 9 

2 

HG06 


BD56604 

22470 t\fo 

3.86 ± 0.04 

150 ± 5 

5 ± 1 

TW 


( 2 ) 

24960 ± 880 

4.33 ± 0.12 

97 

2 

TWHG 



23972 ± 500 

3.999 ± 0.048 

97 ± 10 

2 

HGIO 


HD14321 

21750 ± 130 

3.73 ± 0.04 

150 ± 5 

6 ± 1 

TW 


(3) 

22470 ± 1030 

3.81 ± 0.14 

106 

2 

TWHG 



23601 ± 550 

3.896 ± 0.046 

106 ± 10 

2 

HGIO 


HD14476 

25400 t isg 

3.73 ± 0.04 

92 ± 5 

8 ± 1 

TW 


(4) 

26870 ± 1620 

3.96 ± 0.23 

40 

2 

TWHG 



26006 ± 800 

3.749 ± 0.094 

40 ± 14 

2 

HGIO 


HG2048 

18870 ± 130 

4.27 ± 0.04 

199 ± 5 

4l? 

TW 


(5) 

19530 ±510 

4.30 ± 0.08 

125 

2 

TWHG 



18007± 198 

3.944 ± 0.027 

125 ± 8 

2 

HG06 


HG2255 

19670 ± 170 

3.70 ± 0.04 

329 Ig 

4l? 

TW 


(6) 

21220 ± 760 

3.89 ± 0.12 

318 

2 

TWHG 



20360 ± 178 

3.560 ± 0.022 

318 ± 10 

2 

HG06 


HG2794 

21500 ± 130 

3.77 ± 0.04 

148 +_l 

8 ± 1 

TW 


(7) 

22410 ± 1060 

3.94 ± 0.21 

79 

2 

TWHG 



21624 ± 291 

3.611 ± 0.027 

79 ± 8 

2 

HG06 


W168 

19790 

4.11 ± 0.04 

41 

-13 

1 ± 1 

TW 


(15) 

18870 ± 680 

4.00 ± 0.08 

44 

2 

TWHG 



19259 ± 350 

3.98 ± 0.04 

44 ± 14 

2 

HGIO 


References: TW = This Work; T WHG = This Work u sing H 7 only; 
HGIO = 


Huang, Gies & McSwain 

(20101; and HG06 = 

Huang & Gies 

T^‘ TT' * 

|2006a|b 
















